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We provide a unified graphical calculus for all Gaussian pure states, including graph transformation rules 
for all local and semi-local Gaussian unitary operations, as well as local quadrature measurements. We then 
use this graphical calculus to analyze continuous-variable (CV) cluster states, the essential resource for one- 
way quantum computing with CV systems. Current graphical approaches to CV cluster states are only valid in 
the unphysical limit of infinite squeezing, and the associated graph transformation rules only apply when the 
initial and final states are of this form. Our formalism applies to all Gaussian pure states and subsumes these 
rules in a natural way. In addition, the term "CV graph state" currently has several inequivalent definitions in 
use. Using this formalism we provide a single unifying definition that encompasses all of them. We provide 
many examples of how the formalism may be used in the context of CV cluster states: defining the "closest" CV 
cluster state to a given Gaussian pure state and quantifying the error in the approximation due to finite squeezing; 
analyzing the optimality of certain methods of generating CV cluster states; drawing connections between this 
new graphical formalism and bosonic Hamiltonians with Gaussian ground states, including those useful for CV 
one-way quantum computing; and deriving a graphical measure of bipartite entanglement for certain classes of 
CV cluster states. We mention other possible applications of this formalism and conclude with a brief note on 
fault tolerance in CV one-way quantum computing. 
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I. INTRODUCTION 

The invention of one-way quantum computing (QC) HI 
launched an intensive research effort into this new method of 
QC that eliminates unitary evolution and relies solely on adap- 
tive measurements on a highly entangled state of many qubits 
called a cluster state |2|. Concurrently, other work was un- 
derway generalizing qubit-based QC to QC using continuous- 
variable (CV) systems [3l. The two paths merged with the 
invention of CV cluster states ID, which were shortly there- 
after shown to be capable of serving as the entangled resource 
in the CV version of one-way QC fS'.'Sl. 

Initially, C V cluster states and the platform of one-way QC 
making use of them were not believed to be a promising con- 
tender for scalable QC |5|. It was, however, believed that CV 
cluster states would be convenient for demonstrating the basic 
principles of one-way QC since generating such states in the 
optical context was easier than making ordinary cluster states 
from optical qubits |7, 8|. The main reason for this belief 
was that CV cluster states could be generated deterministi- 
cally, while getting single photons to interact required non- 
deterministic gates whose (heralded) failure happens a large 
fraction of the time |9|. Nevertheless, ideal CV cluster states 
are not achievable since they are singular states (i.e., infinitely 
squeezed) and thus have an infinite average photon number 
and infinite energy. Approximate states must therefore be 
used instead, necessarily leading to errors in any CV one-way 
QC protocol 161 [lOl (we will say more about fault tolerance 
in Section IVl). The most natural choice for these approximate 
states would be multimode squeezed states, but the originally 
proposed method of making them |5| involved experimentally 
arduous inline squeezing operations lfTTHT3]| . This limited the 
expected usefulness of the technology. 

Shortly after the invention of CV one-way QC, it was 
shown that inline squeezing was not required at all and that 



CV cluster states could be generated optically using offline 
squeezing plus interferometry [14]. This method involves 
preparing one single-mode squeezed vacuum state per node 
of the cluster and sending these states through an appropri- 
ately designed network of beamsplitters. In fact, this method 
can be used to make any Gaussian state at all fTsl. This rep- 
resented a vast simplification for experiments, which quickly 
demonstrated the viability of this method of generating CV 
cluster states and their usefulness for simple CV quantum in- 
formation processing tasks 1 16-20|. 

Concurrent with this work was a separate initiative to gener- 
ate optical CV cluster states in a single-shot, top-down fashion 
using just a single optical parametric oscillator (OPO) con- 
sisting of a nonlinear crystal within an optical cavity |21 1. In 
this implementation, independent modes are not separated in 
space (as in previous optical proposals) but are instead taken 
to be the different frequencies within a single beam. The ini- 
tial proposal showed that a single OPO and appropriately de- 
signed multifrequency pump beam could, in principle, gener- 
ate any approximate CV cluster state with a bipartite graphPl 
which includes the universal family of square-lattice graphs. 
Further work revealed that this method could generate a mul- 
titude of small CV cluster states |22 1 or a universal family of 
CV cluster states Il23ll24ll . using a method that has excellent 
scaling potential up to a few thousand optical modes with cur- 
rently available technology I.25-.-.28J . 



' The reader should be aware that we use the term "cluster state" where other 
authors might prefer "graph state," since soinetimes in the literature "clus- 
ter state" is used to refer only to a graph state with a square-lattice graph. 
We would refer to such states as "cluster states with a square-lattice graph." 
Each convention has its proponents, but in the present context, where "CV 
graph state" could have three different meanings (CV cluster state, "H- 
graph state, or general Gaussian pure state), this convention also serves 
a clarifying function. 



Yet another method 1*291 reintroduces the experimentally 
challenging Cz gate. But in this case, only one such gate 
is needed because the modes are encoded temporally, each 
traversing the same optical path (but at different times) and 
each passing multiple times through the same optical hard- 
ware implementing the Cz gate. This method has the addi- 
tional advantage that the cluster state is extended as needed — 
simultaneously with measurements implementing an algo- 
rithm on it — in a manner analogous to repeatedly laying down 
additional track in front of a moving train car (a "Wallace and 
Gromit" approach; see footnote in Ref. |29|). Such a method 
eliminates the need for long-time coherence of a large cluster 
state because only a small piece of the state exists at any given 
time. 

While every CV cluster state — ^regardless of how it is 
made — can be represented by a graph [6J, the single-OPO 
generation method revealed another type of graph that is use- 
ful for describing Gaussian pure states fSOl. This graph in- 
dicates the strength and pairings of the two-mode squeezing 
interactions that act within the OPO, and its adjacency ma- 
trix defines the interaction Hamiltonian directly. Thus, we call 
these graphs Hamiltonian graphs ll2TI424l|3Tl , or H-graphs for 
shortrl Despite this natural way of representing the Hamilto- 
nian interactions by graphs, the resulting states are not "CV 
graph states" in the sense of a CV cluster state with the same 
graph as the "H-graph, although they can be interpreted as CV 
cluster states with (in general) a different graph 1211 . This 
creates an ambiguity in the meaning of a "CV graph state." 

Independently of this work, Zhang showed that ideal CV 
cluster states admit graph transformation rules that correspond 
to local Gaussian operations mapping them to other ideal CV 
cluster states |32|. These rules bear some similarity to the cor- 
responding rules for qubit cluster states 1331 l34ll but they are 
not exactly equivalent. (Related work has also been done for 
qubit stabilizer states and local Clifford transformations ll35l - 
[38l .) Further revisions showed that ideal CV cluster states 
with weighted graphs were necessary for a more complete un- 
derstanding of the graph transformation rules 1 39J — a conse- 
quence of the continuity of the quantum variables in question, 
as opposed to the binary nature of qubits, whose graphs are 
necessarily unweightedP] The effect of quadrature measure- 
ments, which can be used to implement any Gaussian oper- 
ation in CV one-way QC |l6l|4l], has recently been incorpo- 
rated into the formalism, as well B2l . 

This original graphical calculus, while useful for demon- 
strating local Gaussian equivalence of CV cluster states, has 



^ In Ref. 1211 the term "two-mode-squeezing graph" wa.s used instead of '"H- 
graph." These terms are synonymous, and only the latter will be used in 
this paper. 

^ Weighted graph states for qubits have been defined |40 1 but they are not 
stabilizer states. When dealing with CVs, however, weighted graphs occur 
naturally because the entangling operation that makes a CV cluster state 
necessaiily has a strength (which can be — but need not be — chosen equal 
for all interactions). This strength becomes the weight for the con'espond- 
ing edge in the graph. Unlike their similarly named but non-analogous 
qubit counterparts, CV cluster states with weighted graphs are CV stabi- 
lizer states. 



several limitations. First, for the weighted as well as for 
the unweighted case, only ideal (i.e., infinitely squeezed) CV 
cluster states can be represented. As mentioned before, these 
states are not physical. Neither their approximating Gaussian 
states nor any other Gaussian state can be represented in the 
formalism. Second, there are many Gaussian operations (for 
instance, the very common Fourier transform) that do not map 
CV cluster states to other CV cluster states and thus cannot be 
represented as a graph transformation. Third, no connection is 
made with H-graphs; the rules apply only to CV cluster-state 
graphs. Nonetheless, Zhang's formalism is exciting because it 
promises an intuitive visual way of manipulating CV cluster 
states, paralleling similar tools for qubit cluster states ll33l[34l . 

In this paper we generalize these rules in a consistent fash- 
ion to cover all Gaussian pure states, including approximate 
CV cluster states. This includes physical states generated by 
the action of a Hamiltonian with an associated H-graph. Be- 
cause the details of a Gaussian pure state are displayed in its 
graphical representation (and represented uniquely within it), 
this formalism can be used to quantify the deviations from 
ideality for any approximate CV cluster state. Furthermore, 
the formalism can also be used to identify the "closest" CV 
cluster state to any given Gaussian pure state, and it is use- 
ful when considering physical systems whose ground states 
would be useful for CV one-way QC |43 1. We also make con- 
nections with a measure of bipartite entanglement in Gaussian 
pure states. In certain cases, this admits a simple graphical 
rule. 

In what follows, we shall (1) define the unique graph asso- 
ciated with any Gaussian pure state, (2) derive the transforma- 
tion rules for all Gaussian unitary operations, (3) illustrate the 
effect of local Gaussian unitary operations in graphical form, 
showing that they faithfully generalize Zhang's rules |l39l 14211, 
and (4) illustrate the connection of this formalism to approxi- 
mate CV cluster state generation via an H-graph Hamiltonian, 
and (5) provide several applications of the formalism to anal- 
ysis of physical states and the H-graph generation method, as 
well as Hamiltonian ground states and bipartite entanglement. 
The connection with "H -graphs answers an important question 
about the method proposed in Refs. ll2TH24ll — namely, what 
happens when the method is used on physical states. Previous 
connections between H-graphs and CV cluster-state graphs 
have only been rigorously made in the unphysical limit of in- 
finite squeezing. This formalism allows the important effects 
of finite squeezing to be properly accounted for while remain- 
ing entirely within the intuitive framework of the graphical 
representation. 

The mathematics behind this formalism is the complex ma- 
trix formalism for representing and manipulating Gaussian 
pure states [44J. When we interpret these matrices as adja- 
cency matrices for complex-weighted graphs, transformations 
using the symplectic representation can also be interpreted 
in graph-theoretic terms. In this formalism, real-weighted 
graphs, representing idealized, unphysical, infinitely squeezed 
Gaussian states 132] [39l l42l . are generalized to complex- 
weighted graphs that uniquely specify realistic, physical, 
finitely squeezed Gaussian pure states. This extension in- 
cludes generalizing the real-valued graph-state nullifiers for 



ideal CV cluster states f^l to complex-valued nullifiers for 
physical CV cluster states. These results are closely related 
to the stabilizer formalism for Gaussian pure states, which, 
though utilized to some extent for a proof of the CV version 
of the Gottesman-Knill theorem |45i and often mentioned as 
a straightforward generalization from finite-dimensional sys- 
tems 1 46 1, has not been fully explored yet. While we are 
building on existing mathematics, what is new in this work — 
beyond the straightforward graph-theoretic interpretation of 
that mathematics — is showing its natural connection to CV 
cluster states, plus all of the examples, applications, and av- 
enues for future work that open up as a result. We will have 
more to say about the context and importance of our work in 
Section |V] 



II. GAUSSIAN PURE STATES AND 
COMPLEX- WEIGHTED GRAPHS 
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FIG. 1. Original formulation of weighted graphs for ideal CV cluster 
states. Ideal CV cluster states are represented by undirected graphs 
with real- weighted edges f4H6l 1 141 1211124) . (Unweighted graphs are 
a special case with all weights equal to 1.) Each graph uniquely de- 
fines a recipe (i.e., a quantum circuit) for creating a CV cluster state, 
as illustrated above: (1) each node represents a state that is infinitely 
squeezed in the p quadrature |0) ; (2) Cz gates are applied between 
modes in accordance with the graph, with the weight g of an edge 
corresponding to the strength of the interaction Cz{g) = e**''®'' be- 
tween the two nodes connected. These states are unphysical because 
they cannot be normalized. Instead they are approximated in physi- 
cal applications by very highly squeezed states. 



A. CV Cluster States 

The motivation for a graphical study of Gaussian pure states 
begins with CV cluster states L4-^. In the ideal case, CV clus- 
ter states are prepared beginning with a collection of N zero- 
momentum eigenstates, which we write as |0) , where the 
p-subscripted kets satisfy p\s) = s\s) . These states are 
then entangled via a collection of controlled- .Z operations, de- 
noted Cz — sxp{igq (E) q), where the real number g is the 
strength of the interaction. Since all Cz operations commute, 
they can be performed in any order (or simultaneously), which 
leads naturally to the use of graphs as recipes for generating 
particular CV cluster states. An example of such a graph — 
and the CV cluster-state recipe it encodes — is illustrated in 
Figure [T] Each node represents a zero-momentum eigenstate, 
and edges indicate a Cz operation to be performed between 
the two connected nodes. The strength g of the interaction is 
indicated by the label, or weight, of the associated edge. As 
such, weighted graphs with real-valued weights are the natu- 
ral language for depicting ideal CV cluster states. 

Labeling the nodes of the graph in some arbitrary order, we 
can define a symmetric adjacency matrix A = A whose 
(j, fc)* entry Ajk is equal to the weight of the edge linking 
node j to node k (with no edge corresponding to a weight 
of 0). Since a graph is uniquely specified (up to isomorphism) 
by its adjacency matrix, we will often omit the distinction be- 
tween the two and refer unambiguously to "the graph A". The 
collection of controlled-Z operations used to make the CV 
cluster state is now a function of A, denoted C2[A]. The CV 
cluster state associated with the graph A is then 
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where q = (gi, . . . , qn)^ is a column vector of Schrodinger- 
picture position operators. The factor of 2 is necessary be- 
cause each edge weight appears twice in the sum (as Ajk and 
as Akj). 

Ideal CV cluster states satisfy a set of nullifier relations 16] 
I2TI . which can be written concisely as 



(p - Aq) IV^a) - , 



(2.2) 



where p = (pi, . . . , pn)"^ is a column vector of Schrodinger- 
picture momentum operators. This equation actually repre- 
sents N independent equations, one for each component of the 
vector (p — Aq), which are called nullifiers for the state |V^a) 
because that state is a simultaneous zero-eigenstate of them 
(and of any linear combination of them). 

These ideal CV cluster states admit a convenient graphical 
representation in terms of the adjacency matrix A. Some lo- 
cal Gaussian unitary operations |39| and quadrature measure- 
ments Ii42,l can be represented as convenient graphical update 
rules. As discussed in the introduction, this graphical formal- 
ism is elegant and intuitive but also has several limitations — 
most notably, the restriction to ideal (and hence non-physical) 
CV cluster states and to a subset of all local Gaussian uni- 
taries (those which map ideal CV cluster states to other ideal 
CV cluster states). 

Approximate CV cluster states are those for which mea- 
surements of each nullifier give values that are close to zero. 
Quantifying this, we would say that an approximate CV clus- 
ter state is any member of a family of Gaussian pure states, 
indexed by an overall squeezing parameter a for which 



lim cov (p — Aq) = . 



(2.3) 



where the covariance matrix of a vector of operators has com- 
ponents defined as the symmetrized expectation value 
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(2.4) 



with the expectation taken over the state of interest — in this 
case, the a-indexed approximate CV cluster state \iPa{<^))- 
(Notice that this definition assumes that the state has zero 
mean.) The Hmit in Eq. ( |2.3| l is defined component-wise for 
all entries in the covariance matrix, but because the matrix is 
positive definite, it is sufficient to require that the relation hold 
only for the diagonal elements: 

lim (VA(a)l {Pj - V AjkqkY IV-aIc.)) - V j . (2.5) 

fc 

Any a-indexed family of Gaussian pure states {|V'a(«))} sat- 
isfying Eq. p.3| l [or Eq. jl.S) ] defines a family of approximate 
CV cluster states with graph A. We would like these Gaussian 
pure states to be representable directly in a graphical formal- 
ism in a manner that takes into account their unique deviation 
from ideality. 

In addition to this, there are plenty of other useful Gaussian 
pure states that are not approximate CV cluster states at all. 
For instance, a two-mode squeezed state satisfies 



var(gi - (72) = e ^" 
var(pi+p2) =e-^" 



(2.6) 



Such states are readily made in the lab by passing position- 
squeezed and momentum-squeezed beams through a 50:50 
beamsplitter |[3]|47]|48l or directly by nondegenerate paramet- 
ric downconversion ll3l l49l452l . Despite being an a-indexed 
family of multimode squeezed states whose variances tend 
to zero as a — > cx), these states do not satisfy Eq. ( |2.3| l 
for any choice of (finitely-weighted) graph A, which means 
they cannot be represented within the existing graphical for- 
malism for CV cluster states |39| — even in the limiting case 
where a — > cx). This is unfortunate since the two-mode 
squeezed state is related to a two-mode CV cluster state by a 
Fourier transform on one of the modes — a local Gaussian uni- 
tary II2TII that is one of the simplest to perform experimentally. 
That this equivalence cannot be represented in the graphical 
formalism — in either the ideal or the approximate case — is a 
significant drawback. Our formalism addresses all these con- 
cerns. 



B. Desired properties of the grapliical calculus 

In our attempt to generalize the formalism of Refer- 
ences Il39ll42ll . we desire a unified graphical calculus that has 
the following properties: 

1. All Gaussian pure states can be represented uniquely, 
up to phase-space displacements, as graphs. 

2. All Gaussian unitaries — local or otherwise — can be 
represented uniquely, up to phase-space displacements, 
as graph transformations. 

3. All local projective measurements of quadrature opera- 
tors can be represented uniquely, up to phase-space dis- 
placements, as graph transformations. 



4. The representation of a family of approximate CV clus- 
ter states faithfully limits to the standard graph repre- 
sentation of the associated ideal CV cluster state in the 
limit of large squeezing. 

5. The representation of local Gaussian unitaries and pro- 
jective measurements of quadrature operators acting on 
a family of approximate CV cluster states faithfully re- 
produces Zhang's rules 1391 l42l in the limit of large 
squeezing. 

In addition to these requirements, we would also like the 
graphical calculus to be useful for the following purposes: 

• Visualize the entanglement structure of a Gaussian pure 
state. 

• Consider finite squeezing effects within the graphical 
formalism, including their effect on one-way QC using 
an approximate CV cluster state. 

• Make a connection with "H-graphs and their usefulness 
in generating CV cluster states ll2TH23l . including pos- 
sible graph transformation rules between the two types 
of graphs. 

This is not meant to be an exhaustive list. We expect that other 
uses will present themselves as the formalism gets applied to 
actual calculations. 



C. Symplectic representation of Gaussian pure states 

In the present and the following subsection, we shall review 
the complex-matrix formalism B4l for representing Gaussian 
pure states and their transformations among each other, adapt 
it to our notation, and recast it for our purposes. 

It is well known that all A^-mode Gaussian pure states can 
be created by acting on the ground state of N harmonic os- 
cillator with a unitary operation whose Heisenberg action on 
the vector of quadrature operators is a symplectic transfor- 
mation, followed by a phase-space displacement. These are 
sometimes called linear unitary Bogoliubov (LUBO) transfor- 
mations |3|, but we will not use this term. Furthermore, we 
will neglect the phase-space displacement altogether since we 
only desire that our graphical formalism describe the noise 
properties of the state, which do not depend on overall dis- 
placement. 

Stacking q on top of p to form a column vector called x = 

( ? 1 , the Heisenberg action of a Gaussian unitary operation U 
takes the form 



U^iU ^ Si, 



(2.7) 



where S is a symplectic matrix of c-numbers that acts via ma- 
trix multiplication on x as a vector, while If is a unitary opera- 
tor that acts on the individual operators within x. Specifically, 



2N 



U^x.U^Y.^^ 



k-^k • 



(2.8) 



fc=l 



Notice that in general there would be a phase-space displace- 
ment term, which would give x' = Sx + y, but we are ne- 
glecting this. There is a unique S for every Gaussian U, and 
there is a unique U (up to an overall phase) for every sym- 
plectic S. This correspondence can be chosen to faithfully 
preserve composition and map the identity operator I to the 
identity matrix I, thus giving a symplectic representation of 
the Gaussian unitary group |53J. 

The symplectic nature of S is guaranteed because the com- 
mutation relations must be preserved, giving rise to a sym- 
plectic form rj to be preserved by the Heisenberg matrix ac- 
tion. The explicit form of $7 may be deduced by writing out 
the commutation relations for x and requiring them to be un- 
changed under the Gaussian unitary operation. The canonical 
commutation relations [qj,Pk] — iSjk (with fi = 1) can be 
written succinctly as 



I 
I 



in. 



(2.9) 



where the commutator of two operator-valued vectors is de- 
fined as 



T\T 



(§f^) 



(2.10) 



Note that the transpose operation acts only on the entries in 
the matrix (or vector), leaving the actual operators involved 
alonen Requiring that the quadrature-operator commutators 
remain unchanged after the Gaussian operation gives 

iVt^ [x,x'^] 
= [Sx, (Si)T] 
= S[x,x'^]S'^ 
= iSfJS'^. (2.11) 



'' Note on notational conventions. Because we are dealing with operator- 
valued matrices (or vectors in this case), the transpose operation must be 
carefully defined. In fact, we define it in a natural way to simply exchange 
the indices of an operator-valued matrix [(A"'^)j(c = (A)fcj] and leave 
the entries themselves alone. That is, each entry in the matrix — which 
is itself an operator — does not get a transpose applied to it. It is then 
no longer the case that matrix transposition follows the usual distributive 
rule — i.e., (AB)^ ^ B^ A"*" because the operators end up in the wrong 
order. Rather than being a problem, we can use this feature to define the 
commutator-product of two operator-valued column vectors f and s as in 
Eq. ^2.10^ . This has the desired property of forming a matrix whose entries 
are the commutators in question: 



' jk 



■ SkTj 



This is how Eq. j2.9) should be interpreted. Also note that c-number ma- 
trices acting on the vectors within the commutator factor out: 

[Af, (Bs)'^] = Afs'^B'^ - (Bsf'^A'^)'^ 



: Afs^B^ - 

:A[f,sT]B'' 



A(srT)TBi 



The usefulness here stems from the ability to "vectorize" expressions using 
operators in a natural way. This will be useful in what follows. 



Noting that the i's cancel, this equation is exactly the defining 
relation for any 2N x 2N square matrix S to be a symplectic 
matrix with symplectic form $7. The symplectic nature of S is 
therefore guaranteed and required by the need to preserve the 
canonical commutation relations, which themselves play the 
role of the symplectic form O (up to an overall factor of i), as 
shown in Eq. ( |2.9| l. 

A Gaussian pure state (with zero mean) is uniquely speci- 
fied by its covariance matrix. We will write the symmetrized 
covariance matrix for an operator- valued vector as 



'({f^f-}). 



(2.12) 



which accords with Eq. ( |2.4| i if we define the anti-commutator 
product as 



TnT 



{f,s^}:=fs^+(sf^) 



(2.13) 



which is analogous to Eq. ( |2.10| i, and if we require that Her- 
mitian conjugation (indicated by ^) apply only to the operators 
within the vector without transposing the vector itself. We will 
use the notation ^ to indicate transposition of the vector and 
conjugation of its entries: 



(ft)T = (fT)t 



(2.14) 



These caveats are unimportant for our current purposes be- 
cause x = X , but they will be necessary later on when 
we wish to take the covariance matrix associated with 
non-Hermitian operators. Even in those cases, Eqs. \2.A\ 
and ( |ZT2] i still hold. 

Since every iV-mode Gaussian pure state can be obtained 
by acting with a Gaussian unitary operation on the ground 
state of N independent harmonic oscillators, we can use 
the symplectic representation of this operation to parametrize 
these states. To eliminate units in q and p, we will normaUze 
the covariance matrix of the A^-mode ground state to be 



1, 



COVXo 



(2.15) 



where xq is the vector of Heisenberg operators associated with 
this state. This means that vargoj = varpoj = \ for every 
mode j. The Heisenberg operators for any Gaussian pure state 
can be obtained from xq by acting with a symplectic matrix, 
resulting in a covariance matrix of 

covx = cov(Sxo) 
1 



J{(Sxo)t,(Sxor}) 



= S(covio)S'^ 



(2.16) 



(For some thoughts on how to generalize this to mixed Gaus- 
sian states, see Appendix ID]) Since the covariance matrix 
uniquely defines a Gaussian state, by Eq. ( |2.16| l so does SS . 



To be practically useful for our purposes, we need a graph 
representation of S (or more accurately, of SS ) and useful 
transformation rules for this representation. To this end, we 
will decompose S and use the resulting matrix factors to de- 
fine the adjacency matrix for an associated graph. 

There are a number of ways to decompose a symplectic 
matrix, but the one we will be interested in is the following 
unique decomposition for any symplectic S |44|: 



state. Equation ( |2.18| l defines the Heisenberg quadrature vari- 
ables X — S(u ,v)Xo associated with the state in question, 
and Eq. ( |2.19| l gives the (symmetrized) covariance matrix, 
from which the Wigner function, Eq. ( |2.20[ ), may be read- 
ily obtained. If one wishes to work with wave functions, 
then Eq. ( |2.21[ ) can be used. Inversion of these relations to 
find U and V is straightforward. The ground state corre- 
sponds to U = I and V = 0. 



I 
V I 



U-V2 






Ul/2 



X 
Y 



Y 
X 



(2.17) 



where U is symmetric and positive definite (U = U > 0), 
V is symmetric (but not necessarily positive definite), and 
the third matrix is orthogonal and thus irrelevant in the prod- 
uct SS Pj Therefore, while this expansion is unique for a 
given S ll44l . since we only care about SS we can fix X = I 
and Y = and, after multiplying the right hand side above, 
define 



'(U,V) 



U-i/2 

vu-1/2 ui/2 



(2.18) 



Using Eq. ( |2.16[ ), the covariance matrix associated with this 
state is 



.is 






(u,v) - 2 (U''^) (u,v) 



1 /u-i 

2 I VU^i 



U^'V 



U + VU"^V 



(2.19) 



Using this we can immediately write down the Wigner 
function for the state 



1 



X exp 



-N 
= IT GXp 



2^ ^(u,v)^ 



(^(U,V)^j V^(u, 



V) 



(2.20) 



D. Gaussian pure states as undirected graphs with 
complex-weighted edges 



The complex combination U— iV that appears in Eq. p.21[ ) 
is suggestive of a way to unify the two symmetric matrices 
that define a Gaussian pure state. Instead, we will multiply 
this by i and define 



Z:= V + iU 



(2.22) 



for reasons that will become clear shortly. This complex, sym- 
metric matrix Z is only useful, of course, if it has nice trans- 
formation properties under Gaussian unitary operations. In 
fact, this is the case. Defining 



A B 
C D 



(2.23) 



as the symplectic matrix corresponding to the evolution in 
question, if the initial state corresponds to Z as above then 
the new state after the evolution will correspond to B4l 



Z' = (C + DZ)(A + BZ) 



(2.24) 



We will interpret Z as the adjacency matrix for an undirected 
graph with complex-valued edge weights, thus providing our 
graph representation for any Gaussian pure state. A rigorous 
derivation of this relation and of the unique map from Gaus- 
sian pure states to graphs is included in Appendix [A| 



and since the state is pure, we can also write down its position- 
space wave function (up to an arbitrary overall phase) 



V'(u,v)(q) - 7r-^^/4(detU)^/*exp 



^1/4, 



1 



qT(U-*V)q 



(2.21) 



Notice that q, p, and x — ( p ) are c-number column vec- 
tors that correspond to their respective operator-valued coun- 
terparts. Any of these four equations can be used to define 
a Gaussian pure state from any pair of N x N symmetric 
matrices, U and V, with U > ensuring physicality of the 



E. Gaussian graphs from expectation values of observables 

An important operational question is. How can the graph 
for a Gaussian pure state be obtained from the statistics of 
measurements made on the system? To answer it, it is useful 
to consider the covariance matrix from Eq. ( |2.19| l. We can 
immediately extract U from the upper-left block: 



U = (2 COY q) 



1 



(qq 



^T\-l 



(2.25) 



Once we have U, extracting V from the upper-right block is 
straightforward: 



^ Physically, in an optical setting for instance, this term corresponds to a 
passive interferometer, which can be implemented using just beamsplitters 
and phase shifters. These have no effect on the vacuum. 



V = U({q,p^}) 

-2(qq^)-'({q,P^}) 



(2.26) 



Putting these together gives 

Z = V + iU 

^(qq^)-'(({q,P^})+il) 

^(qq^)"'(({q,P^}> + ([q,P^]> 



2 

1 

2< 
(qqT)-l(qpT) 



(2.27) 



Equation ( |2.27[ i shows how to extract Z from the expectation 
values of qjq^. and qjPk, with the latter obtainable from the 
observables {qjPk + Pk^j) using the form on the second line. 



F. Approximate CV cluster states 

The graph representative of a Gaussian pure state defined 
above is, in fact, the natural way to extend the graph repre- 
sentation of ideal CV cluster state to their finitely squeezed 
Gaussian approximants. The canonical method for creating a 
CV cluster state |5 1 is to squeeze all modes as much as pos- 
sible in the momentum quadrature and then to apply Cz[A] 
in accord with a (real-weighted) graph A. An ideal cluster 
state I^a) from Eq. \2.\\ is obtained by taking the limit of 
infinite initial squeezing on all the modes. Let's see what this 
looks like in our formalism. 

The symplectic transformation corresponding to the canon- 
ically generated CV cluster state consists of two parts: the ini- 
tial single-mode squeezing and the controUed-Z operations. If 
we take all modes to be momentum-squeezed such that their 
variance is reduced by a factor of e^"^^ , followed by C2[A], 
this corresponds to a total symplectic transformation of 



I 
A I 











(2.28) 



Comparing this with Eq. p.l8| l, we can immediately read off 
that U — e^^'I and V = A, and we find that 



-2r 



I 



(2.29) 



corresponds to an r-indexed family of approximate CV cluster 
states with graph A since 



lim Zr = A . 

r— foo 



(2.30) 



However, there are many other families of Gaussian pure 
states that fit the bill, including one that will be useful later 
(in SectionlrV]): 

Zq = Jsech2aI + tanh2aA, (2.31) 

which satisfies Xivna-^rx, Zq — A. Figurel2]illustrates the cor- 
respondence. 



G. "H-graph states 

"H-graph states are generated by a multimode OPO pumped 
by a multifrequency pump beam and have a mathematical con- 
nection to CV cluster states, even though the nature of the two 





tanh 2a 
i scch 2q I iscch2o 



2 tanh 2a 




(a) 



(b) 



i scch 2q I i scch 2a 
— tanh 2q 

(c) 



FIG. 2. Complex-weighted graph representation of approximate 
CV cluster states. Nodes no longer specify any particular input state 
on their own and are henceforth colored black to emphasize this dis- 
tinction. Instead, the state represented by the graph is entirely spec- 
ified by the edge weights, which can now be complex. The real part 
of the graph V still has the same interpretation as in the original 
formalism — i.e., a collection of weights for the respective Cz gates 
applied between the linked nodes, while the imaginary part U corre- 
sponds to initial multimode squeezing that only mixes qs and ps sep- 
arately (see text), (a) The graph from Figurefllis reinterpreted in this 
formalism. Because the graph has only real weights, it does not rep- 
resent a physical state (since U ^ 0) — and rightly so since ideal CV 
cluster states are infinitely squeezed, (b) This is the graph for an ap- 
proximate CV cluster state made by the canonical method |5 6|. The 
imaginary-weighted self-loops indicate the amount of initial single- 
mode squeezing applied. In the limit of large squeezing (r — ^ oo), 
these physical Gaussian pure states limit to the ideal graph in (a). 
Notice that in this context, states with large squeezing in p are repre- 
sented by black nodes with a vanishing imaginary self-loop weight, 
rather than by the white nodes of Figure [T] (c) Here is another ap- 
proximate CV cluster state, distinct from that in (b), that also lim- 
its to the ideal case (a) when a — >■ oo. The states represented 
in (b) and (c) are physically distinct, but because they have the same 
large-squeezing behavior, they cannot be distinguished in the origi- 
nal graphical formalism (see Figure [Til. This new formalism allows 
them to be distinguished at the graph level. 



types of graph are different ll2TH24l . In our graphical formal- 
ism, these states correspond to graphs with purely imaginary 
weights. Specifically, 



Z = 



-2qG 



(2.32) 



where the real, symmetric matrix G is the "H -graph for the 
state, and a is a unitless overall squeezing strength. The 
term '"H-graph" refers to the fact that G specifies the (lin- 
earized) Hamiltonian for the OPO that acts on the vacuum of 
the cavity modes in order to generate the state. This Hamilto- 
nian is defined by 

HiG) = i/iK2j Gjfc(a]aJ. - ajCik) 

= ihK{k^Gk^ - a'^Ga) 

= hK{q^Gp + p'^Gq) , (2.33) 

where k is a squeezing parameter per unit time ||2T1 . If this 
Hamiltonian is applied for time t, acting on the vacuum state 
with U = cxp[-|-H(G)<], then a = 2Ki in Eq. ^37\ . 

Most "H-graph states (i.e., those with a full-rank G) cor- 
respond to CV cluster states in the limit of large squeezing 
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if one phase-shifts modes appropriately f2T^. Recent work 
shows that this method can be used, in a scalable fashion, 
to make many small, disconnected square-graph CV cluster 
states from a single OPO |22|. More importantly, a single 
OPO can also be used to make very large QC-universal CV 
cluster states in a scalable fashion Il23ll24l . 

The connection between the "H-graph G and the Gaussian 
graph Z, as indicated by Eq. p.32| i, is through the exponen- 
tial map, which is generally a nontrivial operation on a graph. 
However, when G is self-inverse (G = I), this connection 
simplifies greatly: 



G^ 



I 



Z = i cosh 2al-i sinh 2a G (2.34) 



In this case, the Gaussian graph Z is just a rescaled version 
of G with additional self-loops. We will explore the close 
connection between "H-graphs of this form and CV cluster 
states of the form of Eq. ( |2.31| ) in Section [TV] 



H. Complex nullifiers 

In Reference |6|, the real -valued nullifier formalism is used 
to illustrate the effect of quadrature measurements on ideal 
CV cluster states. (In an optical setting, this corresponds to 
homodyne detection.) Rules for updating the nullifiers under 
such measurements are derived, but they only apply to ideal 
CV cluster states. In addition, the resulting state after the mea- 
surement frequently is only an ideal CV cluster state up to 
phase shifts (if at all) and thus cannot be represented by any 
real-weighted CV cluster-state graph. Zhang has implemented 
these rules as a graph transformation |42|, but the same re- 
striction to ideal CV cluster states still applies. Here we ex- 
tend the nullifier transformation rules from Reference |6| to 
graph transformation rules for quadrature measurements made 
on any Gaussian pure state, including approximate CV cluster 
states, thus generalizing the results from Reference [42]. We 
do this by first generalizing the real-valued nullifier formalism 
to complex-valued nullifiers, which can be used to represent 
all Gaussian pure states. 

The nullifier formalism for CV cluster states, given by 
Eq. ( |2.2| i, can be extended to all Gaussian pure states using 
the simple replacement of the CV cluster-state graph A with 
the Gaussian graph Z: 



(p - zq) ivz; 



= (-Z I) xJ7z |0) 

= Uz{-Z I)Szi|0) 

= ^7zUl/2(p-^q)|0) 
= 0, 



(2.35) 



where we have used Eq. ( |2.18| l to plug in for Sz — S(u,v)> 
and we note that p — iq = — i\/2a, a vector of operators 
annihilating the ground state. Notice that the nullifier vec- 
tor p — Zq is not unique for a given graph since any c-number 
matrix M, acting from the left, will generate a new vector of 
nullifiers that are also satisfied by the state: 



(Mp - MZq) l^z) = . 



(2.36) 



The action of the matrix M represents forming a new nullifier 
vector from linear combinations of the original nullifiers. If, 
in addition, M is invertible, then this new nullifier vector also 
uniquely defines the state. 

We can take the fact that p — iq is proportional to a vector 
of annihilation operations and make the analogy more explicit. 
Let's define 



az := 

a^ - 



V2 



V-'^\p - Zq) , 
.U-V2(p_z*q), 



^- V2' 
which have the usual commutation relations 



(2.37) 
(2.38) 



[az,ag] = iu- 

2 
= 1, 



1/2 



1/2 



Z I 



[x,x 



I 



-z i)in 



I 



U 



U-l/2 
1/2 



(2.39) 



it 



along with [az, a^] = [a^, a^] = 0. In the case of the ground 
state, Z = il, and this expression reduces to the usual result 
(az — a). In the more general case, these operators can be 
used to derive a Hamiltonian for which the associated graph Z 



is the ground state (see Section IV F i 



With this notation in hand, we can calculate 

cov(p - Zq) = cov(-i\/2U^/^az) 

= uV2({4,a^})uV^ 



uV2{(2aT^a-) + ([az,al 
U, 



(2.40) 



where we used Eqs. ( |2.35| l and p.39| ) to obtain the last line. 
For approximate CV cluster states. It's also useful to calculate 
the covariance matrix of just the real part of Z (namely, V): 



cov[p — Vq] = (^V I) (covx) 



-V 
I 



I 



U" 





u 



2 ' 



(2.41) 



where we have used the explicit form for the covariance 
matrix — the second line of Eq. ( |2.19| l. Comparing this ex- 
pression with Eq. ( |2.3[ ), we get a nice interpretation of Z's 
real and imaginary parts: V is the graph of the ideal CV 
cluster state approximated by Z, and |U is the error in the 
approximation — now expressed in quantitative terms as the 
covariance matrix of the nullifiers p — Vq. For this interpre- 
tation to make any sense, of course, U must be small. Since 
U > 0, we can use the trace to say that ^ tr U is the mag- 
nitude of the approximation error in approximating the ideal 
CV cluster state V using Z. This trace corresponds to the sum 
of the variances of each of the nullifiers: 



-trU = y 
2 ^ 

i 



(V'zl (Pj 



fc 



(2.42) 



which should be compared with Eq. p.5| l. 

This extension of the real nuUifier formalism for real- 
weighted, ideal CV cluster states to a complex nullifier for- 
malism for complex-weighted, physical CV cluster states — 
as well as all other Gaussian pure states — corresponds to a 
similar generalization on the level of the stabilizer operators, 
which, in the CV case, describe position and momentum shifts 
in phase space |6, 14 1. As we shall continue focusing on the 
nullifiers here, we give a brief derivation of the corresponding 
generalized stabilizer representation for Gaussian pure states 
in Appendix |C] 



III. TRANSFORMATION RULES 

In this section, we derive the rules for updating a graph af- 
ter a local Gaussian unitary transformation is performed on 
one or more of the modes or when a quadrature measurement 
is made. Since an arbitrary n-local transformation can be ob- 
tained by composing 1- and 2-local transformations, we will 
treat those cases with additional care. We begin by deriving 
the affects of the symplectic transformations on the underly- 
ing matrix Z, and following that, we make connection with 
the Gaussian unitary transformations they represent, and we 
illustrate these transformations graphically. 



A. n-Local Gaussian unitary operations 

Consider a Gaussian pure state on n+p modes. It's graph Z 
can be written in block form as 



Z = 



T R' 
R W 



(3.1) 



where T is n x ti, R is p x n, and W is p x p. The symmet- 
ric matrix T is the adjacency matrix for the induced subgraph 
of Z formed by considering only the n modes in question. 
The other symmetric matrix W represents the induced sub- 
graph corresponding to the untouched modes. R, of course, 
represents the connection between the two sets of nodes. 

Without loss of generality, we can represent an arbitrary n- 
local operation with the symplectic matrix 



^A B ON 

10 

C D 

,0 1/ 



(3.2) 



and 



^ri 



local 



A B 
C D 



(3.3) 



is the symplectic operation for just the n target nodes. Then, 
defining 



J = (A + BT)- 



(3.4) 



we follow Eq. (2.24 1 and write the transformed matrix as 



Z' = 



C 




D 
I 



C + DT DR' 



R 



W 



C + DT DR^ 



R 



W 



A 
I 



BT BR' 





JBR^ 
I 



B 


-1 



(C + DT)J (C + DT) JBR^ + DR^ 



RJ 



W RJBR 



T 



(3.5) 



This matrix must be symmetric, so we know that the upper- 
right block must simply be J R , but let's work it out and see 
why. The first fact we'll need is that (C + DT)(A + BT)"^ 
is symmetric, since it represents a Gaussian operation on a 
valid graph (recall that T is symmetric and that ImT > 0). 
Therefore, 

-(C + DT)JBR^ + DR^ 
= -J'^(C + DT)'^BR'^ + dr"^ 
= J^ [-(C + DT)'^B + J^'^D] R^ 
= J'^[-C'^B-TD'^B + A'^D + TB'^D]R'^. (3.6) 

The properties of a symplectic matrix ll53l include A D — 
C'^B = I and D'^B = B'^D. Therefore, the quantity in 
brackets is equal to I, and the n-local-transformed matrix is 



Z' 



(3.7) 



(C + DT)(A + BT)-i 
R(A + BT)-i 



(A 
W R(A 



BT) trT 
BT) iBR^ 



This has a nice interpretation. The upper-left block is just 
the transformation that would result from applying the n-local 
Gaussian to the nodes in question, without any connection to 
other nodes. The bottom-right block reflects the fact that any 
changes to the other nodes' induced subgraph are mediated 
solely by connections with the target nodes and in a additive 
fashion on W. These are the two key observations — the in- 
duced subgraph on the target nodes cares not about connec- 
tions to other nodes, and the effect on adjacent nodes is only 
through connection with the target nodes. The off-diagonal 
blocks illustrate the action of the transformation on the edges 
connecting the two sets. After discussing the important spe- 
cial cases of 1- and 2-local transformations, we will illustrate 
these graphically to provide additional insight. 



B. Local Gaussian unitary operations 

In the case of (l-)local Gaussian operations, so-called 
LG operations, with 2x2 symplectic matrix 



Slg — 



a b 
c d 



(3.8) 
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this transformation is just 



dt 



t rM ^Slg^ \ a + ht a + bt 

r W / I r bvY^ 

W — 



Z'. (3.9) 



^a 



bt 



a 



bt' 



From these resuhs, we will calculate (below) local graph 
transformation rules corresponding to elementary LG opera- 
tions. 



C. 2-Local Gaussian operations 

Arbitrary Gaussian operations can be constructed out of 
1-local (LG) and 2-local (2LG) operations alone. In fact, 
given the availability of all LGs, only a single fiducial 2LG 
is needed to construct any Gaussian operation |54|. The most 
theoretically simple 2LG is just the Gz gate. This operation 
just adds a real constant — ^proportional to the strength of the 
interaction — to the edge in question. 

Considering the applicability of our results to optical imple- 
mentations, we shall also consider another operation (or actu- 



ally, another class of operations), since the Gz gate is diffi- 
cult to implement experimentally [11]. This class will consist 
of beamsplitter interactions. Specifically, we consider only a 
fiducial type of photon-number-conserving interaction. This 
2LG interaction can be used to model any beamsplitter when 
combined with appropriate phase shifts on the input and out- 
put modes: 



SBs(e) = 



Re 
R. 



fcos 6 — sin 9 

sin cos 9 

cos ( 

sin ^ 




(3.10) 



where sin6' is the amplitude reflectivity of the beamsplit- 
ter [55 1 . This form is particularly simple because B = C = 
and A == D = R^, which gives (A + BT)"^ = RJ = R-9- 
This being the only quantity that affects nodes outside the tar- 
get set, this is particularly convenient. The transformation 
of T is given by 



, _ T _ (cosO -sin0\ (Tu Ti2\ ( cos9 sin9\ 

-M,0±«,e _ l^ sin 61 cose* ) \Ti2 T22) \-s\n9 cos 9 J 

_ /Til cos^ 9 + T22 sin^ 9 - T12 sin 29 T12 cos 29+^ (Tn 
~ i^ T12 cos 261 + i (Tn - T22) sin 261 T22 cos^ 6I + Tn sin' 



sin^ 61 



T22)sin26i 
+ T12 sin 26* 



(3.11) 



The final result for the beamsplitter transformation is: 
T' R«RT' 



Z' = 



Re ( 

: 



w 



(3.12) 



T R' 
R W 



R? 
I 



where T' is defined in Eq. p.ll| i, and we have included the 
second (expanded) form to show the formally simple effect 
that beamsplitter interactions of the form Sbs(^) have on 
Gaussian graphs. (As previously mentioned, modeling a phys- 
ical beamsplitter may additionally involve phase shifts, which 
are local Gaussian operations.) 



D. Quadrature measurements 

According to the rules of Reference ||6l, the first thing to do 
when considering quadrature measurements is to define the 
nullifier that corresponds to the measurement outcome. This 
is the new nullifier that the post-measurement state must sat- 
isfy due to projection onto the measurement basis. (This rule 
applies even if the measurement is destructive.) Next, we 
choose an appropriate invertible matrix M such that the en- 
tries in the new (ideal) nullifier vector Mp — MVq are such 



r 



that only one of them fails to commute with the nullifier corre- 
sponding to the measurement (something which can always be 
done). This new nullifier vector also uniquely defines the pre- 
measurement state (since M is invertible), but because only 
one of its entries (i.e., a single nullifier) fails to commute with 
the measurement nullifier, all the remaining ones will also be 
nullifiers for the post-measurement state as well. The one that 
fails to commute is therefore discarded and replaced by the 
measurement nullifier to form the post-measurement nullifier 
vector. 

Everything said in the previous paragraph remains valid 
when considering the complex nullifiers from Section II H in- 
stead of the ideal ones from Reference [6|. We will there- 
fore focus on one specific example: q-measurements. Such a 
measurement on node j with outcome Sj means that the new 
state has qj — Sj as one of its nullifiers (we assume the de- 
tector noise is negligible). However, since we are neglecting 
displacements, the post-measurement state will instead have 
Qj as the measurement nullifier. The usual nullifier vector for 
a Gaussian graph Z is p — Zq, which is akeady in a fortu- 
nate form since all of the nullifiers commute with qj except 
the jth one. This measurement on an ideal CV cluster state 
corresponds to deletion of the measured node from the clus- 
ter, along with all of its links. The effect is the same on Gaus- 
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sian graphs and can be seen in the nullifier formalism in that 
all references to pj are gone, and linear combinations of the 
post-measurement nullifiers can be used to delete all refer- 
ences to qj, as well (since qj alone is one of the nullifiers in 
the set). As an action on the adjacency matrix Z, this corre- 
sponds to deleting its jth row and column. 

Another useful quadrature measurement is a measurement 
of p. For an ideal CV cluster state, such a measurement deletes 
the corresponding node \m\. preserves the links between neigh- 
boring nodes, up to local phase shifts. Because of the phase 
shifts, the resulting state is often impossible to represent as 
an ideal CV cluster-state graph (although it still has a nulli- 
fier representation), but we can do it easily for approximate 
CV cluster states as a two-step process: (1) perform an in- 
verse Fourier transform on the node to be measured, and then 
(2) perform a q measurement as described above, thereby 
deleting it and its links from the graph. This is not equiva- 
lent in total to a simple disconnection, however, because the 
phase shift generates additional connections in the neighbor- 
hood of the measured node before that node is disconnected. It 
is through this mechanism that a measurement of p preserves 
links in the graph while deleting the measured node |6 1. Mea- 
surements of the general quadrature q cos 6 + p sin 6, which 
can be used to perform Gaussian dyna mic^ on encoded CV 
quantum information in the cluster 161 1411 . are represented 
analogously, with a phase shift by 9 replacing the inverse 
Fourier transform before the q measurement. 



E. Graph transformation rules 

Here we illustrate the transformation rules described above 
as actions on the adjacency matrix Z into rules transforming 
the associated graph for several examples. In all cases, the 
original Z is given in block form according to Eq. ( |3.1[ ), while 
Equations ( |3.9[ ) and ( |3.7| i provide the transformation laws for 
local and 2-local operations, respectively. We will focus on 
representative local transformations, then the 2-local interac- 
tions discussed previously, and finally quadrature measure- 
ments. 

Displacement — The graph rule for displacements in phase 
space is trivial: do nothing to the graph. The Gaussian graph 
only represents the noise properties of the state, which are 
unaffected by overall displacements. 

Local shear — The easiest nontrivial local operation to 
represent in this formalism is a local shear in phase 
space: exp{^gq^). The corresponding local symplectic ma- 
trix is 



9lg 



r(.9) 



1 
9 1 



(3.13) 



resulting in the simple transformation 
2 , Sshe..(g), ft + g V 



W 



(3.14) 



where Z is defined in Eq. ( |3.9[ ). Notice that the only edge 
affected by this transformation is the self-loop on the node 
experiencing the shear (and note that W is not shown): 




t + g 



SahcarCs), 



-)> 




Single-mode squeezing — Squeezing (i.e., reducing the 
variance) in p with squeezing parameter r > is equivalent to 
squeezing in q with squeezing parameter —r and is effected by 
the unitary operator Qxp[— ^r{qp + pq)] — exp[ir(a^^ — a?)] 
and represented by the local symplectic matrix 



»LG 



ozc(r) 



e 
e 



and transforms graphs in a particularly simple way: 

e-''r V^ 



oM^ 



(3.15) 



(3.16) 



Notice that this transformation only affects the edges attached 
to the node being squeezed; there is no effect on the neighbor- 
hood of the affected node (represented by W, not shown): 




e ri 



Ssqucozc{'-) 




Phase shift — The transformations above only affect edges 
attached to the node being acted upon (the "active" node). In 
order to make a change to the W-matrix, however, which rep- 
resents the rest of the graph, we need something else. A sim- 
ple operation that accomplishes this goal and has the advan- 
tage of being very easy to implement in many experiments is 
the phase shift (a.k.a. phase delay), which simply rotates the 
phase plane of the active node by an angle 6. The unitary op- 
erator is cxp[— 16'((7^ + P^)] = e^'®/^ exp(— ifla^a), where 
the overall phase e^*^/^ can be ignored. The corresponding 
local symplectic matrix is just a rotation matrixFl 



'LG 



* phase 



{e)-r- 



COS 6 sin 6 
- sin cos 



(3.17) 



* For universal QC, the ability to measure in a non-Gaussian basis is re- 
quired |6|. In an optical context, this can be achieved through photon 
counting. Because the resulting state is no longer Gaussian, such mea- 
surements are not incorporated into this formalism. 



^ There is some ambiguity in the sign used in the definition of the phase shift. 
The convention we are using is consistent with References Ii21,; '23il24l . 
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The associated transformation is 

sin 9 + t cos 6 



„ Sphasc(0) 

Zi I r 



COS 6 + t sin ( 
r 

cos 9 + t sin ( 



cos 9 + t sin 

sin^rr'^ 

W — 



cos 6 + t sin ( 



(3.18) 



Notice that phase shifts on one node can be used to induce 
additional links within its neighborhood (consider the case 
where 1^12 = in what follows; also note that cg — cos 9 
and sg = sin 9): 




Cfl+tae 



Any single-mode Gaussian unitary operation can be repre- 
sented as a graph transformation by appropriately concatenat- 
ing the rules for squeezing and phase shifting fTS). 

At this point it is useful to mention that all of these rules 
agree with the graph rules derived by Zhang [ 39 1 in the limit of 
Im Z ^ 0, but only when the initial and final graphs remain 
finite in this limit. For example, the Fourier transformFlwhich 
corresponds to the phase shift 



5 phase 



(-f) 



gives 



-r 



-1 

1 



-t-irT 



and is represented as follows: 



(3.19) 



(3.20) 





Wii-t-\j 



12 - i Vir2 



^^22 



W22-t-Vi 



This has no corresponding rule in the ideal limit if t = 0, 
i.e., if the active node has no self-loop, which is the case for 
most ideal CV cluster states of interest. Nevertheless, apply- 
ing the Fourier transform twice has a corresponding rule in the 
ideal limit, even when t — 0, since F~ = —I: 






t -r' 
r W 



(3.21) 



This rule corresponds exactly to Figure 5 in Ref. fW]. Sim- 
ilarly, the rule for local squeezing corresponds to Figure 6 in 
that reference (and is valid for t ~ 0), while Figures 3 and 4 
in that work can also be derived from the rules given here. 

Controlled-Z gate — Similar to the local shear operation 
discussed above, the controlled-Z gate Cz{g) — exp(iggiq2) 
is the easiest 2LG operation to represent in the graphical for- 
malism. The corresponding symplectic matrix is 



*2LG 



.(.9):= 



/I ON 

10 

.9 1 

\.9 h 



(3.22) 



resulting in the simple transformation 



Til Ti2 + g „T 



2 ^ Scz(g)^ lT2i+g T22 
R 



(3.23) 



W 



The graphical representation of this is simply to add g to the 
edge weight between the two active nodes: 



Sc^Ca) 
Wu ' ^ Tia + 





All such gates commute, and thus they can be performed in 
any order or even simultaneously. Despite this theoretical sim- 
plicity, their experimental difficulty 1 1 1 1 suggests developing 
another rule for a canonical 2LG operation. 

Beamsplitter — The beamsplitter interaction Sbs(^) de- 
fined in Eq. p.lO| l, which corresponds to the uni- 
tary exp[-i9{qip2 - hPi)] = ex p[-6>( aiQ^ - a2a{)] and 
whose action on Z is given by Eq. p.ll[ ), has the following 
graph transformation rule (also note that W is unaffected and 
that Cg = cos 9 and sg = sin 9): 



iTll4+T22S^g-Ti2S29) 




SbsCS), 



T12C20 

^(Tll-T22) 
^2 




(^Riicg — Ri2Se) 
Wii 

\Rl2Cg + RilSg) 



{T22ci+Tiisl+Tus2e) 



' Unlike the sign ambiguity for the definition of the phase shift, the Fourier 
transform is fixed by the requirement that measuring in j5 is the same as 
applying the inverse Fourier transform F^ and then measuring in g J5j . 



A useful special case of the above rule is the 50:50 beamsplit- 
ter, for which 9 = j: 
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5(^11 + T22) — Tl2 

-7j(-Rii — R12) 



i(Tii-T22) 




^{Rl2+Rll) 



5(^11 +T22)+Ti2 



Any multimode Gaussian operation can be represented as a 
graph transformation by appropriately concatenating the rules 
for arbitrary beamsplitters, single-mode squeezing, and phase 
shifts through the Bloch-Messiah decomposition 1 15 1. 

Measurement of q, p, and other quadratures — The rule 
for projective measurement of q is to delete the node and its 
links from the graph: 



measure q 
W12 ' ^ 



This is the only measurement rule we need because, as shown 
Section 





m 



HID 



the rule for measuring quadratures other 
than q is to phase shift the node to be measured so that a 
subsequent measurement of q is equivalent to the intended 
quadrature measurement on the original state. For example, 
an inverse Fourier transform F^^ = Sphasc(f ) followed by a 
q measurement implements a p measurement on the original 
state. The graph rule for this is just a concatenation of the rule 
for the phase shift followed by node deletion. 

Because phase shifts can induce new links in the neighbor- 
hood of the shifted node (i.e., it can change W), measure- 
ments other than of q will in general preserve links between 
nodes that were previously mediated by the deleted node. We 
show a special case of this by measuring p on the previous 
graph with W12 = 0; 

W^ii-t-V? 



W^22 



Notice that since the measured node mediated a connection 
between the two other nodes, measuring it in p preserved this 
connection in the form of a new edge connecting those two 
nodes directly. Also notice that the strength of this new con- 
nection is proportional to both of the weights ri 2 in the origi- 
nal mediated connection. These rules agree with Zhang's rules 
for quadrature measurements on ideal CV cluster states 



IV. APPLICATIONS 

A. Closest CV cluster state to a given Gaussian pure state 

Given the many different ways to make approximate CV 
cluster states ll4ll51 [T4ll2T^E4l|291 , a useful question to ask is. 





What is the closest CV cluster state approximated by a given 
Gaussian pure state? Imagine that you are given a system 
with a Gaussian graph Z and you want to know if you can use 



it as a C V cluster state for one-way QC. From Section II H 



we know that V is the graph for the ideal CV cluster state 
approximated by Z and that ^ tr U is the approximation er- 
ror (which vanishes in the ideal case). On the other hand, 
we also know that H-graph states have V = and a diverg- 
ing i tr U in the limit of large squeezing, yet they are useful 
for CV one-way QC l2TH24ll . This naive prescription is there- 
fore not enough. 

In fact, transforming these H-graph states into a useful form 
requires phase shifting nodes appropriately [211, which, of 
course, also transforms the graph Z i-> Z'. The "H-graph 
method, first proposed in Reference |21 1, creates approxima- 
tions to ideal CV cluster states with a bipartite graph. The 
nodes of a bipartite graph can be colored such no edge links 
two nodes of the same color The prescription from Ref- 
erence |21| for using these states requires first performing 
a Fourier transform (i.e., phase shift by — |) on some of 
the nodes and then using the resulting state as an ordinary 
CV cluster state (see Reference [21 1 for details). The trans- 
formed Z' has a nonzero V' and a small approximation error 
(in the limit of large squeezing), which is the basis for the 
main result in that paper. But is this prescription the best we 
can do? Or might there be some other ideal CV cluster state 
obtained by a different phase shifting procedure that is better 
approximated by the original ?^-graph state? 

We will address this question in the next subsection, but to 
answer it, we need to generalize our notion of what it means 
for a Gaussian graph to serve as a CV cluster-state graph. 
Specifically, we must allow the freedom to phase shift any 
mode arbitrarily since this operation can be absorbed into the 
measurement-based protocol to be implemented on the state 
and is thus simply a redefinition of the quadratures and does 
not need to be actively implemented. Our task then is to min- 
imize the approximation error ^ tr U' over all possible phase 
shifts (where the prime indicates the new graph after these 
phase shifts have been applied). Once such a minimum is 
obtained, the resulting V' is the "closest" ideal CV cluster 
state approximated by Z' and thus also the closest one ap- 
proximated by the original Z. 

The details of the calculation are somewhat involved, so we 
relegate them to Appendix |B] choosing to list here just the 
results. Let 



Se 



cos 9 sin 9 
— sin 9 cos 9 



(4.1) 



where 9 = diag(6'i, . . . , ^at). This is the symplectic ma- 
trix representing the phase shifts Oj to be performed on each 
node j. In terms of the original graph Z, the new graph is 

Z' = (-sin6» + cos0Z)(cos6/ + sin0Z)"^ (4.2) 

A sufficient condition for ^ tr U' to be an extremum (see Ap- 
pendix iBli is that 



(ImZ'2),, 



(U'V%--0 Vj 



(4.3) 
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in other words, when corresponding rows of U' and V' are 
orthogonal. A sufficient condition for such an extremum to 
also be a local minimum (see Appendix [B) is that 



Im[(I + Z'2)oZ'] >0, 



(4.4) 



where o represents the Hadamard (entrywise) product of two 
matrices. An interesting generalization of this condition is 
when Im[(I + Z'^) o Z'] > 0, which is the best one can do 
when there is a continuous manifold of phases that all (lo- 
cally) minimize ^ trU'. We will see an example of this in 



Section IV C when we analyze the two-mode squeezed state. 



We would like to stress that these are sufficient conditions 
for local minima only; some minima may not be able to be 
found this way, and not all minima found in this way may be 
global minima. But in certain useful cases (such as the ones 
that follow), we can apply these results to provide evidence 
that we have found the closest CV cluster state (up to phase 
shifts) for a given graph. With rigorous proofs of optimality 
left to future work, our purpose here is to illustrate a useful 
application of this graphical formalism. 



B. Analysis of the H-gjraph method of construction 



As shown in Section [Tl G| all Gaussian pure states created 
using the H-graph method |21 1 necessarily satisfy Z = iU = 
jg-2aG^ where G is the H-graph that defines the multimode 
squeezing Hamiltonian, Eq. ( |2.33| l, that acts on the vacuum, 
and a > is an overall squeezing parameter. To transform 
this state into an approximate CV cluster state, the H-graph 
method prescribes phase shifting particular nodes in accor- 
dance with the desired CV cluster-state graph ||2TI . We will 
thus partition Z into blocks as in Eq. p.l| i in accord with the 
partitioning of nodes to be shifted and nodes to be left alone: 






(4.5) 



We now perform a Fourier transform F on the nodes corre- 
sponding to the upper-left block, using Eq. p.20| i for each or 
applying Eq. p.7|) directly, resulting in 



Z' = 



iu; 



-ur'u^ 



-UaU^^ iUs - iUaU^^Uj 



(4.6) 



We want to know whether this procedure results in a minimum 
in the approximation error ^ tr U' . This is difficult in general, 
but we will do so for a particularly useful case shortly. For the 
moment, we can prove something weaker but still interesting: 
this particular choice of phase shifts results in an extremum 
of i tr U'. This is easily seen since every entry in Z' is either 
purely real or purely imaginary. Therefore (U'V')jj = for 
all j, and thus, by Eq. ( |4.3| l, ^ tr U' is an extremum. 

In fact, this holds for any set of phase shifts by multiples 
of I on purely imaginary graphs, since the above construc- 
tion didn't depend on any particular node(s) being selected for 
phase shifting, and shifting by tt just applies negative signs, 
which doesn't change the real/imaginary nature of the entries. 



Many of these will result in rather large values of the approxi- 
mation error, but multiples of ^ mean that these cases are still 
local extrema. The prescription in Reference |21| explicitly 
chooses the nodes to be phase shifted so that the desired ideal 
CV cluster-state nullifiers vanish in the limit a — > oo. What 
was believed, but not shown explicitly, is that for any finite 
value of a, the prescription minimized these variances. While 
we have still not shown this rigorously in the general case, we 
have provided additional evidence for this claim by showing 
that it results in an extremum of the sum of these variances. 
We will be able to say more about particular examples, pre- 
sented next. 



C. Two-mode squeezed state 

The simplest nontrivial example of "H -graph construction 
of a CV cluster state is embodied by the two-mode squeezed 
state, which results from downconversion in a nondegenerate 
OPO Il3l l49lj52l . This procedure applies the Hamiltonian from 
Eq. ( |233| ) with an "H-graph 121] 



1 

1 



Using Eq. ( |2.32| l, the state created hasj 

Zi 



-2aG 



i cosh 2a 
—i sinh 2a 



~i sinh 2a 
i cosh 2a 



(4.7) 



(4.8) 



where a > is an overall squeezing parameter, and the 
subscript is used because this state will be compared below 
to a canonically-generated two-mode CV cluster state. This 
state is symmetric under exchange of the nodes, so we can 
choose to Fourier transform either one, either of which, using 
Eq. ( |l20l ), results in 



K 



fi sech 2a tanh 2a 
I tanh 2a i sech 2a 



(4.9) 



This is a two-mode CV cluster state with weight tanh 2a, 
which goes to 1 in the limit a — > oo, and approximation 
error ^ tr U'j^ = sech 2a, which vanishes in the same limit. 
This means that 7/^ — > G in this limit, making trivial the con- 
nection between the generated CV cluster state and its gener- 
ating "H-graph for the two-mode squeezed state. This prop- 
erty is not generic — most "H -graphs generate entirely differ- 
ent CV cluster states |21 1 — ^but a particular class of H-graphs 
admit this trivial connection. This is discussed further in Sec- 
tion |IVD] 

For finite a, any combination of possible phase shifts by 
multiples of ^ (see Section IV B i results in an approxima- 



tion error of either cosh 2a or sech 2a. Z'^ corresponds to the 



' Alternatively, the same state can be generated by interfering a g-squeezed 
state Z = ie^" with a p-squeezed state Z = «e~^" at a 50:50 
beamsplitter (3] El ED— i.e., by applying Eq. (3lO\ with g = J to 

7. = fie-^" A 
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latter — a (local) minimum in the approximation error. Let's 
try to use Eq. ( |4.4| i to verify this, though: 

lm[(I + Zf ) o Z'J = 2sech2atanh2 2a K M , (4.10) 

which has eignevalues and 4sech2Q;tanh 2a, and thus 
only the weaker condition, lm[(I + Z'j^) o Z'j^] > 0, is sat- 
isfied. The zero eigenvalue in this case corresponds to the fact 
that any additional phase shift by 9 on the first node can be 
completely canceled by an additional phase shift by —0 on 
the second node. This means that there is a one-parameter 
manifold of phase shifts that all result in Zi i-> Z'j^, which 
has approximation error sccli 2a and which is the minimum 
value obtainable. Specifically, all phase shifts (6*1,^2) sat- 
isfying 61+62 = — f will result in minimal error as an 
approximation to an ideal CV cluster state. In addition, a 
second manifold defined by ^i + ^2 = +f gives a similar 
minimum-error approximate CV cluster state but one which 
has weight — tanh 2a. 

We can also create a version of this state as a canonical CV 
cluster state. In this case, using Eq. ( |2.29| l gives 



1 



le 



-2r 



(4.11) 



Both Z']^ and Z2 are complete graphs on 2 nodes with self- 
loops, and while they are identical in the infinite squeezing 
limit (a — ?► 00, r — >^ cx)), the weightings on each are differ- 
ent for any finite amount of squeezing. While real-weighted 
ideal CV cluster-state graphs fail to illustrate this, the complex 
graphical formalism captures the difference. 

The importance of this difference comes from the relative 
efficiency of each method in its use of squeezing resources. 
It is known that the canonically-generated CV cluster states 
(e.g., Z2) are inefficient in this sense because the resulting 
state has single-mode-squeezed marginals 161 HU. An "effi- 
cient" state has (qjqj) = {Pjpj} and | {{qj,pj}) = for 
all nodes j, which means that all of the correlations are be- 
tween quadratures variables of different systems ll56l . Recall- 
ing Eq. ( |2.20| i for the covariance matrix in terms of Z, these 
requirements become 



(U- 



vu-iv-u-^)„- = (vu-i),, = o 



This fails for Z2 because 

(U2 + V2U2-IV2 - U2-I),, = e-^- , 



Vj. (4.12) 



(4.13) 



which vanishes only in the infinite-squeezing limit (r — ^ cx)). 
In contrast, Eq. ( |4.12[ ) does hold for Z\. Thus, the 'H-graph 
method of constructing a two-mode CV cluster state is effi- 
cient in its use of squeezing resources, while the canonical 
method is not. This fact cannot be seen in the real-weighted 
graphical formalism, but the complex formalism reveals it. 



D. Bipartite, self-inverse ?^-graphs 

Notice that for the two-mode squeezed state created us- 
ing the "H-graph method, Z'l ^- G in the limit a — > 00 of 



large squeezing. This is a general feature of "H -graphs that 
are bipartite and self-inverse [22 1. Such graphs include some 
with square-lattice topology that are useful for universal one- 
way QC |24|. Recalling Eq. ( |2.34[ ), when G is self-inverse 
(i.e., G^ = I), the resulting H-graph Gaussian pure state is 



Z = ie"2"*^ = i cosh 2al - i sinh 2a G . 



(4.14) 



In this case, the Gaussian graph Z is just a rescaled version 
of G with additional self-loops. When G is also bipartite, it 
can be written as 



G = ( ° Go 



where Gg is square and satisfies Gq Gq 
ing 



GnG 



o«^o 



i cosh 2a I 
—i sinh 2a Go 



-i sinh 2a Gq 
i cosh 2a I 



(4.15) 



I, giv- 



(4.16) 



Performing a Fourier transform on the either set of nodes in 
the bipartition 1,21 J . using Eq. p.20| i, gives 



—,_ I i sech 2a I tanh 2a Gq 
l^tanh 2a Gq i sech 2a I 

— i sech 2a I + tanh 2a G , 



(4.17) 



which satisfies Z' — >^ G for large squeezing (a — ?> 00), corre- 
sponding to an ideal CV cluster state with the same graph as 
the H -graph G I22ll24ll . The two-mode squeezed state from 



Section IV C is the simplest special case of this construction. 



We already know from Section IV B that phase shifting by 
multiples of ^ the nodes of a Gaussian pure state created using 
the "H -graph method of construction results in approximate 
C V cluster states with extremal approximation error ^ tr U. 
In this case, the possible values of this error are 



1 



tr U = n sech 2a + (Ni, — n) cosh 2a , 



(4.18) 



where Nb is the total number of nodes in one of the biparti- 
tions (equal to half the total number of nodes in the graph). 
Clearly the minimum of these choices occurs when n — Nb, 
corresponding (nonuniquely) to the original prescription 121]: 
Fourier transform all nodes in one of the bipartitions, and do 
nothing to those in the other set. This results in i tr U = 
Nb sech 2a, which vanishes in the limit a — > oo. 
Let's try to verify Eq. (|4.4|i for this choice: 



Im[(I + Z'2)oZ'] 

= Im [(2 tanh^ 2a I + 2i sech 2a tanh 2a G) 
o [i sech 2a I + tanh 2a G)] 

= 2sech2atanh2 2a( ^ ^^ GJq o Gq 

\Go O Gq I 

= 2sech2atanh2 2a(I + G)o (I + G) > , (4.19) 

where we have used the Schur product theorem (A o B > 
if A, B > 0), and I + G > since G is self-inverse. There- 
fore, Im[(I + Z'^) o Z'] > 0, which we also found inde- 



pendently for the two-mode squeezed state in Section IV C 
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In some cases of interest, like for the two-mode squeezed 
state, there will be a manifold of phase shifts that result in 
equivalent minimal-error CV cluster states made using a bi- 
partite, self-inverse "H-graph, thus accounting for the vanish- 
ing eigenvalue(s) in Eq. ( |4.19[ ). While this does not rigorously 
prove that the extremal value of the approximation error in 
the H-graph construction method for self-inverse, bipartite T-L- 
graphs 1 22, 24 1 is a global minimum (or even, strictly speak- 
ing, a local minimum), these calculations suffice to illustrate 
the usefulness of the complex graphical formalism and sug- 
gest further avenues of research in this area. 



E. GHZ state 

Also discussed in the Uterature lET] [30l El ED is the 
CV Greenberger-Horne-Zeilinger (GHZ) state, which can be 
made for any number of systems N using a complete T-L- 
graplpj with no self-loops: G = I — J, where J is the 
N X N matrix of all ones. We note that the diagonal of G 
is not fixed for this state, but there are restrictions on it. In 
order to obtain a fully squeezed state in the limit a ^- cx), 
an H-graph must be full rank (2T1 (i.e., all eigenvalues must 
be nonzero). In order to obtain GHZ entanglement, self-loops 
on the complete graph must be adjusted so that G has at least 
one eigenvalue of each sign 1.30. .57, .58.1 . Since the spectrum 
of J is {N, 0, . . . , 0), any H-graph will suffice that is of the 
form /3I — J with /? > 0. We choose a hollow G (i.e., zero 
diagonal, corresponding to f3 = 1) for calculational conve- 
nience and because it has been studied specifically in the lit- 
erature Il57ll58l . 

Since J = A^'^'^^J for fc > 1, the Gaussian graph associ- 
ated with G is 



Z^.e-e-^.e-fl+i^^^"^)'^ 



le 



'2a 



N 
e 



fc=i 

2aN _ 2 



fc! 



N 



(4.20) 



As prescribed in Reference ETl . we wish to perform a Fourier 
transform on the first node using Eq. ( |3.20| l. To do this, we 
partition the graph as follows: 



Z = 



(t r r ■ ■ ■ r\ 

r W w ■ ■ ■ w 

r w W ■ ■ : 



w 



(4.21) 



\r w ■ ■ ■ w Wj 



' Alternatively, just like for the two-mode squeezed state (see footnoted, 
the CV GHZ state can be made using single-mode squeezing and interfer- 
ometry [i59J . 



where 



t^W 



N 
i 



-{N + e 



2aN 



1) 



Ne 



N- 



-[l + (iV-l)e 



N-\ 



N 
i 



-ie 



2aN 



1) 



Ne 



N 



-[1 



^A'l 



(4.22) 



(4.23) 



and where we have defined the small parameter e := e~^" 
because we are eventually interested in the limit a — )■ oo. 
After applying Eq. p.20| l, the resulting graph is 

/t' r' r' • • • r' \ 

r' W' w' ■ ■ ■ w 

Z' = r' w' W '■■ : , (4.24) 



\r' w' ■■■ w' W'J 

for which we will evaluate each term exactly and also to the 
lowest nontrivial order in e: 



t' 



W' ^W 



1 



iNe 



N-l 



t l + {N-l)e 



N 



iNe 



Af-l 



-1 



.N 



t 1 + {N - l)e^ 
ie 



r- .2 + {N-2)e 



t 



t 



'1 + {N -l)e 



N 



-1 + iVe^ 
i2e, 



1-e 



N 



'1 + {N ~l)e 



N 



(4.25) 
(4.26) 
(4.27) 
(4.28) 



In the infinite-squeezing limit (a — >^ oo, which corresponds 
to e — >^ 0), all terms vanish except r' , resulting in an ideal CV 
cluster state with a star graph and —1 for all edge weights: 



Z'^ 



/O -1 
' -1 



V-1 



0/ 



(4.29) 



The phase-shifted node is the center of the star A similar 
connection between complete graphs and star graphs is known 
for qubit cluster states Il33l[34l . We hope this formalism will 
be a useful tool in generalizing results like these (for qubits) 
to the realm of CVs. 

That a Gaussian pure state generated from a complete 
H-graph corresponds — after a Fourier transform on one of 
the nodes — to a star-graph CV cluster state in the infinite- 
squeezing limit is already known fT\\, but the existing graph 
transformation rules for ideal CV cluster states |39| do not 
allow representation of the necessary Fourier transform op- 
eration. What is new here — and what this construction 
illustrates — is a unified presentation that includes approxi- 
mate CV cluster states and H-graphs (through the exponen- 
tial map), plus intermediate Gaussian pure states and the rules 
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for transforming between them, which captures all of the de- 
tails associated with finite squeezing wholly from within the 
graphical formalism. 



F. CV cluster states as ground states of a 2-body Hamiltonlan 

Here we derive a (nonunique) Hamiltonian whose ground 
state is a particular CV cluster stata^or, more generally, any 
given Gaussian graph Z. The straightforward Hamiltonian to 
consider is 



H[Z] = 2aiaz 



(4.30) 



which satisfies H[Z] Jt/iz) = by Eq. ( |235| ). This two-body 
Hamiltonian is also positive definite by construction, which 
means that \tpz) is its ground state. Instead, however, we will 
use 

i/[Z] = (p-Zq)H(p-Zq) 

= P^P q^Z*p - p'^Zq + q'^Z*Zq 

= P'^P - q'^(-iU + V)p - p^(iU + V)q 

+ q'^(u2+V2-iUV + iVU)q. (4.31) 

This two-body Hamiltonian, too, satisfies iJ [Z] j-^z) = by 
Eq. p.35[ ), and it is positive definite by construction, which 
means that \ipz) is its ground state. Reference |21| proves 
that any CV cluster state with a bipartite graph is equivalent 
to an "H-graph state up to phase shifts of — ^. Restricting to 
these graphs simplifies the construction even further: 

H\ilJ] = p'^p + iq^Vp - ip'^Uq + q^U^q . (4.32) 

Considering that 

iq'^Up - ip'^Uq = i^iqjUjkPk - PjUjkqk) 
= i^iQjUjkPk -PkUjkqj) 



E 

3 



Ur 



-trU 



(4.33) 



is just an overall shift in energy (which can be ignored), we 
have the result that 



F = p^p + q^U^q 



(4.34) 



' ' After the initial appearance of this work, but before its publication, another 
result appeared addressing this idea 1 60 1 . The analysis is limited to canon- 
ical CV cluster states |5| and uses the original nuUifier formalism |6|, but 
the possible use of ground states of quadratic Hamiltonians for CV one- 
way QC is explored in much greater depth. 



has an "H-graph-representable ground state (with G ex 
— log U). Thus, for any desired CV cluster state with a bi- 
partite graph ETI . some Hamiltonian of this form exists that 
has the desired state as its ground state (up to phase shifts). 



G. Bipartite entanglement 

A general A^-mode Gaussian pure state may exhibit some 
form of multipartite entanglement. In this case, it is useful 
to consider the entanglement for any bipartite splitting of the 
state. Given an A^-mode Gaussian pure state, we split the 
modes into two sets, one with n and the other one with m 
modes, N — n + m. We wish to calculate the entanglement 
between the two sets. 

Since the overall state is pure, we may use the entanglement 
entropy for this, which is simply the entropy of one of the sub- 
systems with the other traced out. For a Gaussian pure state, 
the entanglement entropy is a simple function of the symplec- 
tic eigenvalues ||6T1 of the covariance matrix corresponding to 
the nodes in question. The symplectic eigenvalues are sim- 
ilar to ordinary eigenvalues of a matrix, except that a sym- 
plectic product is taken between the matrix and its symplectic 
eigenvectors instead of the ordinary matrix product and the 
magnitude of the imaginary part is then taken. Thus, while 
the matrix equation Mv^ = XjVj defines the ordinary eigen- 
values of the matrix M, the following equation defines the 
eigenvalues associated with the symplectic product; 



sr2v. 



^jVj 



(4.35) 



Notice the presence of the symplectic form ft in this rela- 
tion. Further note that the matrix SJl is not Hermitian, since 
we have (Sn)^ = fl'^T, = -rJS, but both SO and fJS 
have the same eigenvalues (since they are related by a simi- 
larity transformation). Therefore, for an A^-mode matrix, we 
obtain 2N imaginary eigenvalues \j± — ±iaj, which oc- 
cur in conjugate pairs. The N symplectic eigenvalues are just 
the a-j > 0. 

We want the symplectic eigenvalues of the reduced density 
matrix, though, so we have to consider the covariance matrix 
truncated to the subset of n modes. Let P be a matrix con- 
structed in the following way: first, create a diagonal matrix 
that has a 1 in the diagonal entries corresponding to the nodes 
of the set to be kept and O's everywhere else; then, remove the 
all-zero rows. Also, let P — ( g p ), which is also not a square 
matrix. Here, P is an n x A^ matrix, and so P is a 2n x 2N 
matrix. 

The covariance matrix of the reduced state of n modes is 

- - T 

then S := PSP , a 2?! x 2?! matrix. Keep in mind that 
this truncated covariance matrix does not have a correspond- 
ing graph in our formalism, because it is generally not pure 
(if it were pure, then the original state would be trivially a 

- - T 

product state). Since O = POP is the symplectic form for 
the truncated set of modes, we seek the ordinary eigenvalues 
ofPOSP^. 

Our goal is to find a simple way to use the graph Z to 
read off the bipartite entanglement entropy across an arbitrary 
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boundary dividing the graph into the two subsets of nodes. In 
particular, we would like to find an interpretation of the sym- 
plectic eigenvalues, and hence the entanglement in terms of 
the shape and the weights of the given graph. For a general 
graph and arbitrary divisions into subsets of n and m modes, 
this interpretation is not so straightforward. However, for spe- 
cific graphs and bipartite splittings, a simple, instructive con- 
nection between the entanglement and the graph can some- 
times be found. 

A particularly straightforward example is that of the canon- 
ical CV cluster states fS^. Canonical CV cluster states are of 
the form Z = iel + V, where e = e"^*". As shown in Fig- 
ure [2] the imaginary self- loops represent the initial squeezing 
of each node, while the matrix V contains the real weights be- 
tween neighboring nodes (corresponding to the strength of the 
Cz gates, which also contain squeezing 1151 ). The covariance 
matrix for such a state is 



S := 



1 /I 



V 



2e \V eH + V 



(4.36) 



This means that in order to obtain the ordinary eigenvalues of 

- - ^ 

PJISP , we need to solve 



= det 



V-AI eH + W' 
I V AI 



det(A2l - V% e^j _^ ^) ^ 



In other words, besides the initial squeezing r, the entangle- 
ment between any node k and the rest of the graph is deter- 
mined by the number of its neighbors and the strength of the 
link with each neighbor. Only those "nodes" of the graph that 
have no neighbors at all would give a ~ ^, corresponding to a 
pure reduced state and hence no entanglement with the actual 
nodes of the graph. For any given r, however, any link with 
nonzero weight means cr > ^ and thus entanglement. 

Quantitatively, both an increasing number of neighboring 
nodes and an increasing strength of the links enhances the 
entanglement, because the entropy grows with a, and (y ^ \ 
represents the mean excitation number of the reduced thermal 
state. For the special case of equal unit weights Vki — 1, we 
have (T — \\/l + e^^'^Nk, where Nk represents the number 
of neighbors of node k. Thus, the maximum entanglement 
obtainable between a single node and the rest of the graph 
is determined by the maximal number of neighboring nodes, 
the so-called connectivity C :— maxfc Nk- For constant con- 
nectivity, like in a 2D lattice with C = 4, the entropy will 
be bounded and does not increase with the size of the lattice. 
Only for an increasing connectivity do we get larger entropies, 
and in principle the entanglement between a single node and 
the rest of the graph may grow unboundedly with the number 
of its links. This result is consistent with the bounds on the 
offline squeezing per node needed to create a canonical CV 
cluster state of constant connectivity (61. 



(4.37) 



where V = PVP"^ and W = PV^P"^, and the second 
line follows because the bottom two blocks commute. No- 
tice that since I — P P is a projector, W — V = PV(I — 
P P)VP > 0, and thus we will label the eigenvalues 
of W - V by t-l > 0. Equation ( |437| i then gives A| + e^ = 
— fc'? The 2n eigenvalues iij± of flT, are therefore given by 



(4.38) 




which gives us the n symplectic eigenvalues of S: 
1 




(4.39) 



,~,2 



The eigenvalues i^^ of W — V contain the information about 
the graph in question. In particular, when we consider the 
bipartite entanglement between any single node k of the graph 
and the remaining nodes (n — 1, m ~ N — 1), we have 



v^v^Vkk, w^w = Y,v^^ 



kl ■ 



(4.40) 



As aresult, we simply have to solve A^+e^ ~ — J2\i\i^k} ^ki' 
leading to a single symplectic eigenvalue 



V. CONCLUSION 

We have generalized weighted graphs for continuous- 
variable (CV) cluster states in a natural way to a graphical 
calculus for all Gaussian pure states. The mathematics behind 
this generalization is not new |44, 45, 53 62 1. What is new 
is interpreting the matrix formalism of Simon, Sudarshan, and 
Mukunda |44| as transformations on an undirected, complex- 
weighted graph. This would be a mathematic triviality if it 
were not for the fact that applying this graphical interpreta- 
tion to approximate CV cluster states and letting the overall 
squeezing go to infinity results in exactly the same graphs as 
are already being used to represent ideal CV cluster states (4\- 
I6l[l4, 21-241. In addition, the graph transformation rules im- 
plied by the formalism immediately generalize all of the ex- 
isting graph transformation rules ll32l l39l l42l to any Gaussian 
pure state and limit to these rules in the (unphysical) case of 
ideal CV cluster states. It is these remarkable facts that make 
these results important. 

This graphical formalism satisfies all five of the essential 



desired properties outlined in Section II B We have also made 
headway on the three bullet points that followed. Our achieve- 
ments with the formalism thus far include using it to 

• incorporate all details of finite squeezing within the CV 
cluster-state graph (Section[lIF[); 



--Jji+e+*^E^.v 



{i\i^k} 



(4.41) 



• distinguish between different approximants to a given 
ideal CV cluster state at the graphical level (Sec- 
tion irTFli; 
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incorporate "H-graphs ll2TH24ll within the same graph- 
ical formalism through the exponential map (Sec- 
tion HTgIi; 



• generalize the nullifier formalism ||6l to all Gaussian 
pure states (Section [lIH| l; 

• use the nullifier formalism to define both an error matrix 
and a scalar approximation error for an approximate C V 
cluster state (Section [llH| i; 

• define matrix and graphical transformation rules for a 
complete set of simple Gaussian unitary operations and 
quadrature measurements (Section[lIl|; 

• define the "closest" ideal CV cluster state to a given 
Gaussian pure state (Section IV A|i; 



analyze the optimality of the H-graph construction 



method with respect to this notion (Section IV B i, 
including the specific examples of the two-mode 
squeezed state (Section IV Q and a useful subclass 
of H-graphs (Section|IVD|l; 



• demonstrate generation of a star-graph approximate C V 
cluster state from an approximate GHZ state made us- 



ing an H-graph (Section IV E i; 



• identify classes of two-body Hamiltonians that have CV 
cluster states as their ground state (Section [IV F| i; 

• quantify bipartite entanglement in terms of the graph- 
ical formalism, with the explicit example of canonical 
CV cluster states (Section|IVG|. 

We anticipate that this list will grow over time. Specifically, 
we expect that the formalism will serve well when consider- 
ing the propagation and manipulation of quantum information 
through approximate CV cluster states using homodyne detec- 
tion and development of computer software for visualizing the 
effects of Gaussian operations on CV cluster states, as well as 
other uses not yet discovered. 

We conclude with a few words about the prospects of fault 
tolerant quantum computing using Gaussian approximate CV 
cluster states. Recent work |10| has demonstrated that all 
Gaussian approximants to ideal CV cluster states are inher- 
ently faulty when used for one-way quantum computing sim- 
ply due to the fact that they are finitely squeezed. This has 
led to suggestions (made in private communication) that these 
results forbid a fault-tolerant implementation of this quantum 
computing paradigm. This is not the case. 

That C V cluster states are error-prone has been known since 
the beginning [SJ . The main conclusion that should be drawn 
from Reference ifTOll is that there is no "magic pill" to elim- 
inate the effects of finite squeezing using a simple qubit (or 
other) encoding scheme. Instead, fault tolerance must be ad- 
dressed from the very beginning because unlike qubit clus- 
ter states, which remain physical in the limit of zero errors 
in preparation and use, CV cluster states are unphysical in 
this limit, since ideal states require infinite energy. This is 
not in any way a show-stopper, however, and the authors of 



Reference 1 10] go to great lengths to show several possible 
approaches to error correction that do not fall victim to their 
no-go theorems. (It would be interesting to apply the matrix- 
product-state calculations from that paper to noisy qubit clus- 
ter states to see if analogous results are obtained in that con- 
text, in order to ensure we are comparing apples with apples.) 
Fault tolerance in CV one-way quantum computing remains 
an important open problem and an active area of research. 
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Appendix A: Derivation of the complex- weigh ted graph 
representation 

Here we derive the connection between a Gaussian pure 
state and its complex-weighted graph Z. There is a vast math- 
ematical literature on the set of all allowable graphs Z, known 
as the Siegel upper half-space. The main results for our pur- 
poses can be found in Reference | |53]| . with more details in 
Reference ll44l . An extensive review of the Siegel upper half- 
space and its connection to the symplectic group can be found 
in the Ph.D. thesis by Freitas ll62ll and the references therein. 

Definition. The graph corresponding to a Gaussian pure 
state is Z :— (qq )~"'^(qp ). 



This graph was defined in Eq. ( |2.27[ ). We wish to prove 
several properties of all such graphs. 

Theorem A.l. Every graph corresponding to a Gaussian 
pure state is complex-weighted, undirected, and unique and 
has positive definite imaginary part. 

Proof. Let Z be a graph corresponding to a Gaussian pure 
state. That Z is complex weighted and unique is guaranteed 
by the definition: expectation values of operator- valued matri- 
ces result in matrices of complex numbers, corresponding to a 
complex-weighted graph. It is unique for a given state because 
expectation values are uniquely determined by the state. 

An undirected graph has a symmetric adjacency matrix. For 
showing the symmetry of Z, the second line of Eq. ( |2.27| i is 
useful, and we repeat it here for reference: 



Z = ^(qq 



T\-l 



*I+({q,pT}) 



(Al) 



To prove the symmetry of Z we refer to Eq. ( |2.16| ), which 
shows that any Gaussian pure state has a covariance matrix 
that is i times a symplectic matrix: 



2S = SS'^ e Sp(2iV, 



(A2) 
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since it is the product of a symplectic matrix S and its trans- 
pose. While S is not unique for a given Gaussian pure state, 
S is, and we partition it as follows: 



S = 



1 f A B 
B^ D 



(A3) 



Being a covariance matrix, S = S > 0. (Being the co- 
variance matrix for a valid quantum state requires more than 
this |3|, but we don't need it.) This implies immediately that 
A = A'^ > 0, and D = D'^ > 0. That 2S is symplectic 
requires, additionally, that 



I 
-I 



J7 = 4I]f2S 



A B 



B^ D 



I 
I 



A B 
B^ D 



AB' BA AD B^ 

JB^)2-DA B^D DB 

From this we can see immediately that 

(BA)'r==BA, 
(DB)T==DB, 
AD = I + B^ 



(A4) 



(A5a) 
(A5b) 
(A5c) 



From the definition of the covariance matrix, we have that 
(qq^) = A, and |({q,p^}) = B. That A = A"^ > 
guarantees immediately that Im Z — A exists and is also 
symmetric and positive definite. The only remaining item to 
prove is that ReZ = A^^B is symmetric. Equation ( |A5a| i 
gives 

(BA)^ = BA 
=^ A"^(BA)'^A"^ = A^^BAA^^ 

=^ (A^^B)^ = A^^B, (A6) 



since A = A . Therefore, Z is an undirected graph. 



D 



Theorem A.2. Every complex-weighted, undirected graph 
with positive-definite imaginary part represents a unique 
Gaussian pure state (up to an arbitrary overall phase). 

Proof. Let Z be the graph in question. We'll split it up into its 
real and imaginary parts as usual: Z = iU + V, where U — 
U > 0, and V == V .To find the Gaussian pure state that 
Z represents, we will construct a covariance matrix from it and 
prove that it satisfies the conditions to be a valid covariance 
matrix for a Gaussian pure state. 
We now define 



^^ 1 / u-i u-^ 

^ '^ 2 \Y\]-^ U + VU^^V 



(A7) 



just as in Eq. p.l9| l. In this case, though, we need to prove that 
it is a valid covariance matrix for a Gaussian pure state given 
the assumptions made about Z. A straightforward way to do 
this is to show that it is the result of conjugation of the ground 



state covariance matrix ^I by a symplectic matrix (which nec- 
essarily represents a Gaussian unitary operation). We define 



Sz:= 






(A8) 



paralleling Eq. ( |2.18| l. One can verify directly that Sz^JSz = 
Jl, guaranteeing the symplectic nature of Sz. Then we have 



;szs 



z^z : 



(A9) 



which shows that Z represents a valid Gaussian pure state. 

To show uniqueness of the state represented (up to overall 
phase), we assume that there exists another covariance ma- 
trix Sz 7^ Sz that is represented by Z. Recalling Eq. ( |A5c| i, 
which holds for an arbitrary pure-Gaussian-state covariance 
matrix, we write the blocks of Sz as 



Sz = 



1 



B^ A 1 



B 

f A^^B^ 



(AlO) 



Using Eq. ( |A1| |, the graph representing this state is A^ (zl + 
B). By assumption, this must be the same as Z = iU + V, 
which gives A = U^^ and B = U^^V. Plugging these back 
into Eq. ( |A10| l shows that Sz = Sz, in contradiction with the 
assumption. Therefore, Sz is unique for a given graph Z. D 

Theorem A.3. The transformation law for graphs represent- 
ing Gaussian pure states under Gaussian unitary operations 
is given by Eq. ( |2.24[ ). 

Proof. Since we have a unique way of passing from the co- 
variance matrix for a Gaussian pure state to its graph and back 
again, our method will be to show the action of an arbitrary 
symplectic operation on the covariance matrix and then ex- 
tract the new graph from it. Rather than dealing with the co- 
variance matrix alone, it will be useful instead to consider the 
combination 



n 



-i ""' 
2 \y\]- 

1 f\ 



-1 



U"^V-iI 



2 VZ 

2 ^Z 

I 
Z 



U" 



-1 fi zn 



1 



2i 
I Z 



(Z-Z 



vu^^v 



(I z* 

-i -1 



(I z* 



(All) 



and similai-ly for Sz' — f^i, with Z i-> Z'. Equation p.24[ ) 
requires that Sz' — SSzS , with 



S = 



A B 
C D 



(A12) 
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from Eq. ( |2.23| l, repeated here for reference. Using the sym- 
plectic property of S allows us to write 



■^Z' 



n 



S Sz--f2 S^ 



= iS 



(I z*)^fl 



Notice that 
I 



I Z*)S^flS 



A + BZ\ 
C + DZJ 



I Z* 
-1 



(I Zns'^. (A13) 



(A + BZ) , (A14) 



where Z = (C + DZ)(A + BZ)-i. The fact that ImZ > 
guarantees that the inverse exists (unless B = 0, in which case 
invertibility of S guarantees that A^^ exists). Equation p.24| i 
then amounts to the claim that Z' = Z and also that Z satis- 
fies all the requirements for representing a Gaussian pure state 
(i.e., symmetry and positive-definite imaginary part). 

To show the symmetry of Z, we use a trick similar to that 
used in Eq. ( |A13| i: 

= Z Z^ 



I Z^U2 



= d z"^) s'^os 



-T, 



= (A + BZ)'^(Z-Z )(A + BZ) 



(A15) 



;:.T 



Since (A + BZ) is invertible, Z = Z . Similarly, to show 
that Im Z > 0, we calculate 



Ini Z 



1 



t(Z 



= 1(1 z*)sTnsQ 



= — (A + BZ)"(Z - Z*)(A + BZ) 
2i 

= (A + BZ)"(ImZ)(A + BZ). (A16) 

Inverting this relation verifies that (Im Z > 0) =^ (Im Z > 

0)- 

We have shown that Z satisfies the requirements for repre- 
senting a Gaussian pure state. All that's left to show is that 
Z' = Z. To do this, we first take the conjugate transpose of 
Eq. ( |A14[ ), giving 



(I Z*) sT = (A^ + Z*BT) (^i z*) 



(A17) 



Plugging Eqs. ( |A14[ ) and ( |A17| i into Eq. ( |A13| l and canceling 
the appropriate factors gives 



^--^^-u 



I Z O 



I z 



(A18) 



Clearly, Z appears everywhere that Z' should appear. The 
reader can check that solving this equation for Sz' and ex- 
tracting its graph Z' = (qq )~^(qp ) does, in fact, show 
that Z' = Z. Since the graph for a Guassian state is unique, 
this verifies Eq. p.24| ). This transformation law is called a 
generalized Mobius transformation, and the interested reader 
is directed to Reference li62J for a more in-depth mathematical 
analysis. D 



Appendix B: Derivation of the closest CV cluster state to a given 
Gaussian pure state 

There are other matrix models for Gaussian pure states be- 
sides the one we are using here. One of these, based on the 
Siegel disc [62], is useful for these calculations. (We will 
forgo presentation of the entire model, referring the interested 
reader to Reference [62J and the references therein.) Based on 
this model, we define 



K := {I + iZ){l-iZ) 



(Izl) 



I -il) 



(Bl) 



We will also need the following: 

I + K = [(I - iZ) + (I + iZ)]{I - iZ)-i 
= 2(I-*Z)-i, 

as well as 



(B2) 



(B3) 



Notice that since I + K is symmetric, so is K. We repeat 
Eq. (|4.1|l here for reference: 



S« = 



cos 6 sin 6 
- sin cos 



(B4) 



where = diag{Oi, . . . ,9^)- We can derive the trans- 
formation law for K with respect to these phase shifts 
[Cf. Eq. ( |Al4l )]: 



K' = (I ^I) f „\ 



I -*I 



(I tl) Se I 2 
-w ,. ..^ /'I 



I *I 



(I -^l)^e[^ 
' (I -*I) 



-ie 



Ke 



'ie 



(B5) 
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An extremum of ^ tr U' occurs when dj tr U' = 0, 
where 9, := tJ-. Let's calculate the left-hand side: 



= ldjtv{l-iZ')+c.c. 

= 9jtr[(I + K')"^]+c.c. 

= tr{-(I + K')-^[dj{l + K')](I + K')-^} + c.c. 

= tr{-(I + K')~^(-i)e-*^^Ke-'^(I + K')"^} 

+ trjtranspose} + c.c. 
= 2i tr{(I + K')"^e-^^^Ke-'^(I + K')"^} + c.c. 

(B6) 



2i[K'{l + K')-\j+c.c. 



where "c.c." stands for "complex conjugate," and all entries 
of 6j are zero except for the (j, j)*, which equals Oj. We will 
also need the following: 

K'(I + K')"^ = (I + K')~Ml - (I + K')"^] 
= (I + K')-^-(I + K')-' 
^1 .„„ 1 

2* 

1 



-(I-^Z')--(I-^Z') 



l\2 



Then 



Ki- 



5,trU'^-(I + Z'2),, 



(B7) 



c.c. 



= -Im(Z'2 



-(u'v; 



03 
jj ■ 



(B8) 



Setting this to for an extremum verifies Eq. ( |4.3[ ). 

To show that we have a minimum rather than just an ex- 
tremum, the Hessian of tr U' must be positive definite. The 
Hessian matrix has entries 



dkdjtrV = -dkiI + Z'^)jj+c.c. 

- '-dk[{I + iZ'){I - i7l)\,, + c.c. (B9) 

In addition to E q. (|B3| l, we can find another similar relation 
by inverting Eq. (|Bl[): 



I + K 



and thus, also. 



= [(I + iZ) + (I-iZ)](I + zZ)-i 

= 2(I + iZ)-S (BIO) 



I + iZ = 2fI + K 



-i\-i 



(BID 



and similarly for the primed matrices. We now define the 
placeholder matrix 



Q':=(I + K'-i)(I + K') 



21 + K' 



K'. 



(B12) 



Notice that 

I + Z'2 =4Q'"^ 
The partial derivatives of Q' are given by 



(B13) 



5fcQ' = 6)feK' + 9feK'-i 

= 9fcK'-K'-i(afeK')K'-i 

= -i7rfe(K' - K'-i) - i(K' - K'-i)7rfe , (B14) 

where the (fc, /c)* entry of TTfc equals 1, while all others are 0. 
We can plug these results into Eq. ( |B9| i: 

dkdj tr U' 
= 2itr[7rjafeQ'"^] + c.c. 
= -2itr[7rjQ'-i(9feQ')Q'"^] +C.C. 
= -4tr[7rjQ'-Vfc(K' -K'-i)Q'-i] +C.C. 
= -4tr{7r,Q'-Vfe[(I + K')Q'-' 



(I + K'-^)Q'-^]}+c.c. 



4tr{7r,Q'-Vfe[(I + K 



i-\\-\ 



(I + K')-^]}+c.c. 



-2 trlTTjQ'" Vfe[(I + ill) - (I - iZ')]} + c.c. 



= -4«tr(7rjQ' V^Z') + c.c. 
= 2Imtr[7rj-(I + Z'2)7rfcZ'] 
= 2Im[(I + Z'2),,Z',^.] 

= 2Im[(I + Z'2)oZV, 



(B15) 



where o represents the Hadamard (entrywise) product of two 
matrices. Requiring the matrix with these entries to be posi- 
tive definite verifies Eq. \\.A\. 



Appendix C: Stabilizers for Gaussian pure states 

To find the stabilizer operators for the finitely squeezed, 
canonical CV cluster states, we start by constructing the sta- 
bilizer of the vacuum state |0) of a qumode hk |45|. For 
the dimensionless quadrature operators q and p, where a = 
4=(g + i'p), we obtain 



|0) = exp(aafe) |0) = exp ^(4 + ipfc) |0) 



(CI) 



Further, from this we need the stabilizer for a single-mode 
squeezed state Sirk) |0), with a squeezing parameter r/j > 
and a ^-squeezing operator S'(rfc) = (i^^\\rk{qkVk +Pfe'?fe)]- 
The stabilizer equation can be written as 

Sivk) |0) = 5(rfe)exp [^(q,, + zpfe)] s\Tk)S(rk) |0) 

= exp [^(e+'-^gfe + ^e-'^-pfc)] S{rk) |0) . (C2) 
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In the case of momentum squeezing, with S{~rk), we have This result corresponds to the complex nullifier 



S{-rk)\0)=S{-rk)exp 
X S{-rk) 



V2 



(qk+ipk) S\-rk) 



= exp 

Let us rewrite this as 
exp 



W)Xk 



M'-"''i^ 






-le+'^'Pk] 



■V2' 



Si-rk) |0) . 
(C3) 



(C4) 



formally using the WH shift operators X{s) = e ^'^^ and 
Z{s) — e*'*''. Now we define a := —^/2e~^^s such that the 
momentum-squeezed stabilizer becomes 



exp I 



ie-2'-^s2 



)Xk{s)Zk{ie--'^''s). 



(C5) 



In the limit of infinite p-squeezing rj, — ?► oo, this operator 
approaches Xk{s), which stabilizes the zero-eigenstate |0) , 



withXfe(s)lO), 



Pfc 



|0) for all s e M, as expected. 



Now we can proceed to create CV cluster states in the 
canonical way: by pairwise applying the Gz gates, indicated 
as C|? for a link between nodes k and I. The N stabiliz- 
ers of the initial N momentum-squeezed modes, showing in 
Eq. ( |C5| l, with /c = 1, 2, . . . , iV, are then transformed for each 
interaction with neighbor / as 



exp ( — |e 



i„-2..^2)C|'Xfe(s)Cfc|'Zfe(ie-2'--s)Cf 



= exp (-ie-2-^s2) Xk{s)Zi{s)Zk{ie-^'"'s) . (C6) 

Eventually, collecting all these interactions, we obtain the N 
new stabilizers 

Xk{s)Zk{ie-^'-'^s)\{Zi{s), (C7) 



exp I 



~2rfc„2- 



leN'ik) 



where JV{k) is the set of neighbors of fc. In the limit of infinite 
squeezing (r^ — > oo), we get back the well-known, ideal CV 
cluster-state stabilizers for unweighted graphs. However, this 
time, the above stabilizers also do the job for finite squeezing 
and uniquely represent the corresponding approximate clus- 
ter state. The nuUifiers are obtained by taking the log of the 
stabilizers: 



exp(-ie-2'-'=s2 
= exp( 


)Xkis)Zk{ie~^'^'s)Y[Ziis) 
iej\f(k) 

— ^e^^''''s^) exp —is (^pk — ie 


-'"■'qk)] 


X exp(+ie"^'''=s^) ||exp(isg/) 

ieM{k) 


= exp 


-isipk - ie ^'■''qfc ~ )> ]qi) 
I 


, (C8) 


for all fc = 1, 2, . . . , iV and for all s e M. The nullifiers are 
therefore 


Pk 


-ie '^'^''qk -)> ]qi Vfc. 
ieM{k) 


(C9) 



(p - Zq) l^z) - , 



(CIO) 



as expressed in Eq. ( |2.35| l, with a complex adjacency matrix 
Z having imaginary diagonal entries ie~'^^'' and the remain- 
ing entries being either or 1 depending on the particular CV 
cluster state with unweighted edges. For example, for two 
modes this reproduces the result in Eq. ( |4.11| i for a canonical 
two-mode CV cluster state. More generally, the result corre- 
sponds to complex- weighted graphs including self-loops. 

Any A^-mode Gaussian pure state can be built from A^ 
squeezed vacua sent through passive linear optics (modulo 
phase-space displacements) lITSl . In terms of stabilizers, this 
means that, without loss of generality, the stabilizers of N 
momentum-squeezed states are transformed as 

exp (-ie-2'''=s2) UXk{s)U^UZk{ie-^'''' s)U^ 

^exp[-is{pi-ie-^-^q[)] , (Cll) 

where p'^, and g^ are the linearly transformed momentum and 
position operators, respectively, after the corresponding (in- 
verse) unitary transformation U. Provided this U represents a 
Gaussian transformation, the exponent on the right-hand-side 
of Eq. ( |C11[ ) will always be linear combinations of the gener- 
ators q and p. This would include the canonical Cz interac- 
tions, as discussed before. However, now we shall restrict our- 
selves to only passive, number-preserving unitaries U, with- 
out loss of generality [ 15 1. The canonical case would then re- 
quire that the squeezing parts of the Cz gates be absorbed into 
the offline momentum squeezers, corresponding to Bloch- 
Messiah reduction |T4l|. 

For the case of a passive linear transformation, we can write 
a = La, where L := X + iY is an N x N unitary ma- 
trix (with real X and Y), and therefore [Cf. the last matrix in 
Eq. ^AT\ ], 



X Y 
Y X 



(C12) 



Finally, through Eq. ( |C11[ ), we arrive at the new stabilizers 

exp|-is^[(Yfc, -ie-^''>'Xki)qi 

+ (Xki + te-^'-'-YkOpi]} . (C13) 



For the nullifiers, we obtain 



where 



(Ap-Bq)|V'z)=0, 

Aki = Xki + ie-^'^'Yu , 
-Bki = Yki - le-'^^'-Xki , 



(C14) 



(C15) 
(C16) 



which we may rewrite as 
(A^^Ap - A-iRq) l^z) = (P - Zq) |i^z> = , (C17) 
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with Z := A~^B. This gives us again the complex adja- 
cency matrix for an arbitrary pure Gaussian TV-mode state. 
We note that there are at most 4iV^ parameters to determine 
the stabiHzer/nulHfier ||45l . These, however, are not indepen- 
dent, as L must be a unitary matrix, and B follows from A. 
A general Gaussian unitary transformation has 2N'^ + N free 
parameters, without displacements, which is the same number 
for representing a symplectic transformation from Sp(2iV, M). 
For representing pure Gaussian TV-mode states (modulo dis- 
placements), it is enough to apply a general Gaussian uni- 
tary transformation to an A^-mode vacuum state, where after 
Bloch-Messiah reduction IfTSJI . the first passive transformation 
has no effect on the vacuum fT4^. Thus, N real squeezing 
parameters Vk and iV^ parameters for the remaining passive 
transformation L suffice to uniquely determine the matrices A 
and B, and hence the state through Z. 

As an example we refer to the standard two-mode squeezed 
state, obtainable by interfering a ^-squeezed state with a p- 
squeezed state at a 50:50 beamsplitter ||3] |42l |48j . Then we 
have 



A 
B 



1 

\/2 V«e 
1 

71 



1 

-2r2 



1 



(C18) 
(C19) 



and 



A-' = 



J_ (1 

x/2 



3+2r2 
+2r2 



(C20) 



resulting in the matrix Zi from Eq. ( |4.8| l using equal initial 
squeezing ri = r2 = a. Another simple example is the A^- 
mode vacuum state, for which rj. = 0, A = L, and B = zL, 
so that Z = il for any L. The vacuum always remains an 
uncorrelated graph with only self-loops. 



Appendix D: Mixed Gaussian states 

There is a very simple special case of mixed Gaussian states 
for which the entire pure-state graph calculus presented in this 
article follows through as well almost trivially. This special 
case is usually referred to as the iV-dimensional isotropic os- 
cillator 1441 . The covariance matrix for a general mixed Gaus- 
sian TV-mode state is given by 



H'il.^-- 



(Dl) 



generalizing the expressions from Eq. ( |2.16| l, by replacing the 
A''-mode vacuum/ground state 1/2 by 1/2 times the above di- 
agonal matrix with K = diag(K), where the vector k = 
(ki, ..., ktv)""" contains the symplectic eigenvalues (times 2) 
of S. For pure states, we have K = I. 

Now let us assume that all symplectic eigenvalues are equal 
to K, corresponding to an A^-mode Gaussian state built from 
N thermal states with identical excitation number, replacing 
the N initial vacua in Eq. ( |2.16[ ). Carrying along this one extra 
parameter k, we can use our graphical formalism to describe 
this special case, as well. We leave to future work a possible 
extension of our formalism to general Gaussian mixed states. 
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